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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  (metric)  units  as  follows: 


Multiply 

By 

To  Obtain 

cubic  feet 

0.02831685 

cubic  metres 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  kelvins* 

feet 

0.3048 

metres 

square  feet 

0.092290304 

square  metres 

*  To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings,  use  the  following 
formula:  C  =  (5/9)  (F  -  32).  To  obtain  kelvin  (K)  readings,  use:  K  =  (5/9)  (F  -  32)  +  273.15. 


AN  INFORMAL  MONOGRAPH  ON 
RIVERINE  SAND  DUNES 


I.  INTRODUCTION 

The  formation  and  behavior  of  sediment  waves  produced  by  moving  fluids 
(principally  water  and  air)  are,  in  equal  measure,  intellectually  intriguing;  of  great 
engineering  importance;  and  aesthetically  engaging.  Because  of  the  central  role  they  play  in 
river  hydraulics,  fluvial  ripples  and  dunes  have  received  extensive  analytical  attention  from 
engineers  for  at  least  the  past  two  centuries  (e.g.,  Sainjon's  1871  publication,  cited  by 
Leliavsky  (1955);  and  Du  Buat's  1786  work,  cited  by  Graf  (1971));  and  even  more 
intensive  descriptive  study  for  even  longer  from  geologists  (formerly,  before  this  age  of 
specialization,  known  as  natural  philosophers).  A  comprehensive  review  and  assessment 
of  the  resulting  body  of  literature  would  be  far  beyond  the  scope  of  this  investigation. 
Indeed,  such  an  effort  would  be  a  scholarly  enterprise  requiring  many  years-perhaps  a 
whole  career— of  effort  by  one  versed  in  subjects  ranging  from  descriptive  geomophology 
to  mathematical  theory  of  autocorrelation  analysis.  The  review  reported  here  was  limited  to 
the  analytical  aspects  of  riverine  bed  forms,  principally  large  dunes  and  bars— those  of 
sufficient  size  to  interfere  with  navigation  and  other  uses  of  rivers.  Herein  these  are 
referred  to  collectively  by  the  geological  term  "megadunes."  Attention  is  focused  on  the 
mechanics  of  megadune  formation,  and  the  geometrical  and  kinematical  characteristics  of 
the  resulting  bed  forms. 

The  modem,  and  far  more  successful,  era  of  bed-form  analysis  was  initiated, 
without  question,  by  the  work  of  Exner  (1925),  which  brought  (unknowingly)  small- 
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perturbation  stability  theory  to  bear  on  the  problem.  A  second  rebirth  of  the  subject  was 
brought  about  by  Nordin  and  Algert's  (1966)  application  of  time-series  analysis  (also 
referred  to  as  autocorrelation  analysis,  or  spectral  analysis)  to  the  description  and  analysis 
of  fully  developed  bed  forms.  Accordingly,  this  study  was  limited  to  Exner's  (principally 
because  of  its  historical  importance)  and  subsequent  developments;  i.e.,  to  the  period  since 
1925. 

The  principal  objective  of  the  study  was  to  ascertain  if  it  is  now  possible  to  make 
reliable  predictions  of  the  conditions  (combination  of  flow,  fluid,  and  sediment  properties) 
under  which  megadunes  will  form;  and  of  the  lengths  and  heights  of  the  resulting  bed 
forms.  By  way  of  preview,  the  answers  turn  out,  unfortunately,  to  fall  somewhere 
between  "no"  and  "very  approximately." 

Chapter  n  of  this  repon  summarizes  and  critiques  the  principal  publications  on  the 
subject,  divided  into  three  categories:  analytical  (principally  stability-theory)  models; 
empirical  relations  resulting  from  dimensional  and  other  types  of  analysis  of  field  and 
laboratory  data;  and  statistical  (principally  time-series  or  spectral)  models.  Chapter  III  is 
given  over  to  a  discussion  of  the  current  status  of  the  problem,  and  summarizes  what  is 
known  and  reviews  the  major  stumbling  blocks  to  further  progress.  In  Chapter  IV  a  new 
equation,  based  on  a  somewhat  unconventional  theory,  for  dune  height  is  developed. 
Verification  of  five  of  the  leading  empirical  dune-height  relations  is  presented  in  Chapter  V. 
Chapter  VI  recommends  a  new  approach  to  dune  management,  and  the  summary  and 
conclusions  are  presented  in  Chapter  VII. 

The  following  considerations,  which  are  somewhat  philosophical  in  nature,  may  be 
helpful  to  the  understanding  of  the  analyses  reviewed  below,  and  to  appreciation  of  the 
difficulties  encountered  in  analysis  of  sedimentary  bed  forms.  First,  it  should  be  borne  in 


mind  that  bed  forms  are  the  result  of  an  instability  that  is  inherent  in  the  flow,  or  more  pre¬ 
cisely  stated,  in  the  interaction  of  the  flow  with  the  bed.  There  appear  to  be  at  least  two 
types  of  instability.  The  first,  which  occure  at  relatively  low  Froude  numbers,  is  an 
instability  that  originates  at  the  flow-bed  interface,  and  is  present  also  in  the  case  of  very 
deep  flows  and  flows  without  free  surfaces  (e.g.,  sediment-transporting  flows  in  closed 
conduits).  It  produces  ripples,  dunes,  and  bars.*  The  second  principal  type  of  instability 
results  from  an  interaction  between  standing  waves  on  the  water  surface  and  the  bed,  and 
produces  anddunes.  The  latter  type  of  instability  will  not  be  considered  further  herein. 

The  first  instability-the  one  that  produces  ripples,  dunes  and  bars-is,  in  certain 
limited  respects,  similar  to  that  which  produces  turbulence  in  fluid  flows.  Each  occurs  only 
over  a  limited  range  of  conditions.  For  other  conditions— sufficiently  low  or  high  sediment- 
transpon  rates  for  bed  forms;  sufficiently  low  Reynolds  number  for  turbulence-the  flows 
arc  stable,  resulting  in  a  flat  sediment  bed  and  laminar  flow,  respectively.  Both  turbulence 
and  bed  forms  can  be  prevented  for  flows  in  the  ranges  of  instability,  provided  almost 
heroic  measures  are  taken  to  prevent  the  flow  and/or  bed  surface  from  becoming  even 
slightly  disturbed.  Any  disturbance  rapidly  amplifies,  and  equilibrium  turbulence  or  bed 
forms  soon  emerge.  Bed  forms  and  turbulence  both  are  comprised  of  a  range  of  fre¬ 
quencies,  and  there  is  flux  of  wave  content  among  the  components  of  different  frequencies. 
In  the  case  of  turbulence,  the  flux  is  from  the  lower  frequencies  (big  eddies),  at  which  the 
turbulence  is  generated,  to  the  higher  frequencies,  where  the  turbulence  energy  is 
dissipated.  Bed  forms,  on  the  other  hand,  arc  generated  at  tht  higher  frequencies  (shorter 
wavelengths),  which  merge  (due  to  the  frequency  dependence  of  celerity)  to  form  waves  of 
lower  frequency  which  are  dissipated  by  the  flow  (Jain  and  Kennedy  1974).  Both  bed 
forms  and  turbulence  in  their  fully  developed  state  are  characterized  by  equilibrium 

*  The  nomenclature  utilized  herein  is  that  adopted  by  the  ASCE  Task  Committee  on  Bed  Forms  in  Alluvial 
Channels  (1966). 
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distributions  of  frequencies.  For  turbulence,  the  equilibrium  spectral  content  varies  as  f*^^ 
(f  •=  frequency),  over  the  equilibrium  range  of  the  well  known  Kolmogoroff  spectrum. 
Equilibrium-range  bed-form  spectra  vary  as  f'^  (Jain  and  Kennedy  1974;  see  also  papers 
cited  therein  for  further  description  of  frequency  content  of  bed  forms). 

Finally,  and  perhaps  most  important,  is  the  fact  that  both  turbulence  and  bed  forms 
are  an  inherent  characteristic  or  property  of  most  flows  of  engineering  interest  in  their 
respective  fields.  Either  can  be  suppressed,  but  only  over  short  distances  or  times,  and 
only  by  extraordinary  measures  (e.g.,  turbulence  suppression  by  flow  constrictions  like 
those  found  upstream  of  wind-tunnel  throats).  Both  have  decidedly  good  effects,  and 
without  them  Earth  would  be  a  far  different  place.  Without  turbulence,  it  could  not,  for 
example,  support  life  as  we  know  it.  Without  bed  forms,  river  stages  would  vary  much 
more  widely  with  rising  and  failing  discharge,  thereby  exacerbating  flooding  and 
navigation  problems. 

Both  turbulence  and  bed  forms  present  engineering  with  formidable  problems. 
Because  the  underlying  instabilities  cannot  be  avoided,  both  turbulence  and  bed  forms,  and 
the  untoward  effects  they  produce,  must  be  allowed  for  and  managed  to  the  extent  possible. 
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II.  SUMMARIES  OF  RESULTS 


A.  Introductory  Remarks.  This  review  presents  the  significant  results  of  the 
principal  studies  in  summary  form.  The  studies  are  divided  into  three  categories  according 
to  the  approach  followed:  analytical,  empirical,  or  statistical.  The  smdies  are  presented  in 
chronological  order  in  each  category.  Each  summary  utilizes  the  mathematical  symbols 
and,  as  appropriate,  some  of  the  sketches  appearing  in  the  source  reference.  The  source- 
reference  figure  numbers  and  figure  captions  are  retained. 

B.  Analytical  Models.  Eleven  studies  of  this  type  were  reviewed,  starting  with  the 
classical  study  of  Exner  (1925). 
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B.l.  Exner  (1925)  (source  reference;  Leliavsky  (1955)). 


Basis  Qf  ihgQK 

1 .  "Erosion”  equation  (sediment-condnuity  and  sediment-transport  relations). 


0v 

dt  dx 


(1) 


2.  Fluid  continuity. 


b(z-q)  V  =  Q  =  const. 


(2) 


3.  Fluid  momentum 


3v  dv  ,  1 

—  +  v  —  +  kv-y  +  — ^  =  0 
dt  dx  p  dx 


(3) 


Principal  results  and  observadons 


1 .  Combining  (1)  and  (2)  yields 


M  dq 
dt  (2-11)2  dx 


dz 

after  neglecdng  — .  General  solution  of  (4)  is 
dx 


(4) 
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(5) 


z-'n=f[-^-x] 

(2-T1)2 


For  initial  profile  Tio  at  t  =  0 


Tlo  =  Aq  +  Ai  cos  a  X 


(6) 


bed  profile  at  time  t  is 


Ivft 

T]  =  Ao  +  Ai  cos  a[x  -  - — -] 

(z-ti)2 


(7) 


These  profiles  are  shown  for  different  t  in  Fig.  B.  1 .2.  Their  celerity  is 


c  = 


M 


(z-T1)2 


(8) 


and  their  amplitude  remains  constant. 

2.  Including  effects  of  fluid  dynamics,  by  combining  (3)  with  (1)  and  (2)  yields, 
after  linearization. 


at2  'dxdt  at  ax2 


(9) 


For  Tio  given  by  (6),  bed  profile  at  any  t  is 


T]=Ao  +  Aie‘^ '^^‘cosalx- 


(10) 
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If  friction  is  neglected  (k  =  0)  in  (10) 


Tl  =  Ao  +  Ai  cos  Ct(x  -  (Ve  +  N2/4  t].  (11) 

"Velocity  of  the  bank  movement"  given  by  (10)  is 

c=|(|-i)  (12) 

c  increases  with  k  (i.e.,  with  decreasing  wavelength,  X). 

3 .  If  friction  is  included,  wave  amplitude  decreases  as  waves  move  downstream. 
For  zero  friction,  amplimde  remains  constant  but  waves  distort  as  they  move 
(see,  e.g.,  figure  B.1.2). 

4.  Because  formulation  is  quasi  one-dimensional,  analysis  is  limited  to  very  long 
bed  forms  (i.e.,  dunes  and  bars;  not  ripples). 

5 .  No  relations  for  wavelength  or  height  of  bed  forms  result  from  this  analysis. 

6.  Analysis  is  noteworthy  as  first  of  stability  type  analyses,  and  for  laying 
groundwork  for  most  future  analyses  of  this  genre. 

7 .  Exner  extended  analysis  to  include  effects  of  variation  of  channel  width. 
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Figyres 


1^ 


Fic.  12.  Key-Uiupram  for  Profeaaor  Exncr'a  nntiition. 


Figure 
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Fxo.  in.  Professor  Exxier  s  mathexsaticaily  derived  dime  prodle. 


Figure  B.1.2.  Profiles  at  different  times  given  by  (9)  and  (10). 
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Notation 


Some  notation  is  defined  in  the  figures. 


AoAl  = 

b 

E 

f 

g 

j 

k 

M 

N 

P 

Q 

t  = 

V  = 


Vm 


defined  in  (6) 
channel  width 

eg 

any  function 
gravity  acceleration 

k 

function  of  a  and  k;  0  <  2  -  j  <  depending  on  X. 
coefficient  in  linearized  friction  term 
b 

2 

gO/bv^j  •  Vm 

pressure 

discharge 

time 

local  mean  (depth-averaged)  velocity 
average  velocity  over  whole  flow 


a  =  Ik/K 

y  =  component  of  fluid  weight  along  channel 
e  =  "erosion  coefficient" 

X  =  wavelength 

p  =  fluid  density 
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B.2.  Anderson  (1953). 


Basis  of  theory 


1 .  Sediment  continuity. 


3Q'  ^  « 
-^  +  — =  0 
9x  9t 


(1) 


2.  Sediment-transport  relation. 


Q'  =  l^u(x,0)  (2) 


3.  Potential-flow  description  of  stationary  water-surface  waves  of  sinusoidal 
form  on  flow  over  a  sinusoidal  bed  (see  Fig.  B.2.1) 


4.  Karman's  logarithmic  velocity  distribution  is  used  to  estimate  u  near  the 
bed,  which  is  required  to  obtain  Q*  from  (2): 

u  =  ^U  (3) 


Principal  results  and  observations 

1 .  Bed  profile  is  given  by 

^  =  cos^mh 

where 
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(5) 


P  Uni  cosh  mh 

~  m  “  2  sinh  mh 

2.  Wavelength  relation  is  obtained  by  equating  amplitude  given  by  (4)  with  sin 
Pt  =  1  to  that  for  equilibrium  flow  over  a  sinusoidal  bed: 

rnhdanh  mh  -  ^  (6) 

Verification  of  (4)  is  presented  in  Fig.  B.2.2. 

3.  Bed-form  celerity  is  obtained  by  substituting  (3)  into  (2),  and  then  the 
resulting  expression  for  7k  into  (5): 

P  Cm  coth  mh  ^ 

c  =  —  = - — —  u  (7) 

™  2YgBVi 

4.  Observations: 

(a)  Anderson  did  not  allow  for  movement  of  bed  forms  in  formulation  of 
the  velocity  potential.  His  formulation  actually  is  for  moving  bed 
waves  passing  under  stationary  water  waves.  This  gives  rise  to  sin  Pt 
term  in  (4). 

(b)  Wavelength  relation  is  obtained  by  equating  maximum  amplitude 
reached  by  bed  form  migrating  under  stationary  surface  waves  to  bed 
amplitude  which  produces  stationary  surface  waves  of  fixed  amplitude 
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a.  Setting  sin  Pt  =  1  in  (4)  is  not  warranted,  and  good  agreement 
shown  in  Fig.  B.2.2.  is  considered  almost  fortuitous. 

(c)  Form  of  celerity  relation  is  generally  correct,  with  c  increasing  with 
decreasing  X  and  increasing  G. 
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Figores 


Fig.  U  -  Definition  Sketch 
for  Evaluating  Length  of  Sediment  Tfavea 

Hgan  B.2.1. 


Froudi  Numb«r  ' 


* 


Fig.  5  ^  Relative  Wave  Length 
as  a  Function  of  Froude  Nianber 


Figure  B.2.2 


Notation 


B 

c 

C 

Fr 

G 

m 

O' 


u(x,y) 

U 

P 

Yg 

7k 

X 


Some  notation  is  defined  in  the  figures. 

=  8.5  +  5.75  log  y/dg 

=  bed-form  celerity 

=  Chezy  coefficient 

=  sand  size 

=  Froude  number 

=  total  volumetric  sediment  discharge  per  unit  width 

=  2it/k 

=  volumetric  discharge,  per  unit  width,  of  sediment  panicipating  in  bed-form 
migration 
=  time 

=  hcKizontal  velocity  at  (x,y) 

=  mean  velocity 

=  cm  =  temporal  frequency  of  bed  fcnrms 
=  bulk  specific  weight  of  bed  sediment 

=  constant,  equivalent  to  Exner’s  "erosion  coefficient" 

=  wavelength 
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B.3.  Kennedy  (1963) 


Basis  of  theory 

1.  Potential-flow  formulation  of  free-surface  flow  over  a  sinusoidal  bed 
moving  downstream  with  celerity  Ub  (see  Fig.  B.3,1.). 

2 .  Sediment  continuity 


dx  3t 


3 .  Power-law  sediment  transpon  relation,  with  phase  shift  5 


G(x,t)  =  m^: 


dx 


(x-5,  -  d,t) 


(1) 


(2) 


4.  Assumption  that  observed  bed  forms  arc  those  with  wavelength  that 
produces  fastest  amplitude  growth;  i.e.,  those  for  which 


^2a 

dtdk 


(3) 


Principal  result  and  observations 

1 .  Wavelength  relation  is  (see  Fig.  B.3.2) 

p2  _y^  1  kd  tanh  kd  +  k5  cot  kS 

8*^  (kd)2  +  (2  +  k5  cot  k5)  kd  tanh  kd 


I 


(4) 
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2. 


Bed-fonn  celerity  for  wavelength  given  by  (4)  is 


T  T  nGk  (■  sinh  2kd  -t-  2kd  1  , 

**  2B  [sinh^kd  -  jkd  cotjkd  -  ij 


(5) 


3 .  Analysis  yields  mathematical  basis  for  classification  of  bed  forms  (see  Fig. 
B.3.3). 


4.  Observations: 

(a)  Phase  shift  5  between  local  sediment  discharge  and  local  near-bed 
velocity  is  of  central  imponance  to  the  theory.  Such  a  phase  shift 
unquestionably  exists,  but  it  has  not  been  well  formulated.  In  actuality, 
this  phase  shift  includes:  phase  shift  between  local  bed  elevation  and 
local  near-bed  velocity;  phase  shift  between  local  bed  elevation  and 
local  bed  shear  stress;  "adjustment  distance"  of  local  sediment-transpon 
rate  to  local  flow  conditions;  and  gravitational  effect  of  local  bed  slope 
on  local  sediment  discharge. 

(b)  This  is  a  linearized  (small  wave  amplitude)  theory,  and  therefore  is 
likely  of  limited  value  in  analyzing  fully  developed  bed  forms. 
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Shon-crested  antidunes 


Figure  B.3.3. 


Notation 


Some  notation  is  defined  in  the  figures. 

a  =  bed-form  amplitude 

B  =  bulk  specific  weight  of  bed  sediment 

d  =  mean  flow  depth 

F  =  Froude  number 

G  =  sediment  discharge  per  unit  width,  on  a  weight  basis 

G  =  average  (along  channel)  value  of  G 

j  =  8/d 

k  =  2rt/L 

L  =  wavelength 

n,m  =  constants  in  sediment-discharge  relation 

t  =  time 

Ub  =  bed-form  celerity 

x,y  =  space  coordinates 

5  =  phase  shift  between  local  sediment  discharge  and  local  near-bed  velocity 

Ti(x,t)  =  bed  elevation 

<(i  =  velocity  potential 
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B.4.  Kennedy  (1969) 


This  review  article  presented  further  development  and  interpretation  of  Kennedy’s 
(1963)  theoretical  model.  The  principal  elements  of  Kennedy's  theory  are  depicted  in  Fig. 
B.4.1.  Corrections  to  the  theory  made  by  Reynolds  (1965)  were  included,  and  the  resulting 
revisions  to  the  predicted  conditions  for  occurrence  of  different  bed  forms  were  made.  The 
principal  new  development,  oudined  below,  was  a  predictor  for  heights  of  bed  forms. 

Basis  of  theory 


Modification  of  Kennedy’s  1963  theory  gives  bed-form  celerity  as 


TT  /m  .  U  1  -  F2kd  tanh  kd 
Ub  =  (Tb/P)nk  u  ,  _  p2kd 


(1) 


2 .  Average  sediment  discharge  in  bed-form  migration  is 

L 

Tb  =  (Ub/L)ao  J  (1  +  sin  kx)dx  =  UbUo  (2) 

0 

Principal  results  and  observations 

1 .  Revised  diagram  for  occurrence  of  different  bed  forms,  given  in  Fig.  B.4.2. 

2.  Limiting  wavelength  relations  were  modified  in  accordance  with  Reynolds’ 
(1965)  findings,  with  results  shown  in  Fig.  B.4.3.  Kennedy’s  basic 
wavelength  relation  ((4)  of  Section  B.3)  remained  unchanged. 
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3 .  (1)  and  (2)  yield  the  following  relation  for  ripple  steepness 


2ao  P  U  -  Uc  tanh  kd  -  F^kd 
^  nji  U  1  -  F^kd  tanh  kd 


Good  agreement  was  obtained  between  the  bed-wave-steepness  relation  and 
flume  data  for  flows  over  coarse  beds  (see  Table  B.4.1). 
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0-8-L/4  .  0*jkd«-r/2 
Antidunes  Moving  Upstream 

(a) 


L/4-  8.L/2,  r/2«jkd«» 
Antidunes  Moving  Downstream 

(b) 


Deposition 


Scour 


h  8 

uU,-d,f) 

L/2<8*3L/4,  T<jkd<3»/2 
Ripples  or  Dunes  Moving  Upstream 

(not  observed  to  occur ) 

(c) 


3L/4<8«L,  3*’/2<jkd<2ir 
Ripples  or  Dunes  Moving  Downstream 

(d) 


Fig.  2.  Mechanisms  of  instability  for  principal  bed  forms  predicted  by  the  poten¬ 
tial-flow  stability  analysis  (u^^,  the  x  component  of  the  velocity). 


Figure  B.4.1.  Graphical  depiction  of  key  elements  of  Kennedy's  (1963)  theory. 
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i 

Fig.  4.  Conditions  for  occurrence  of  various  bed  forms. 


Hguie  B.4.2. 


BorMn  a  Lin 
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0440  mm  Sond 

Plofo 

0949  mm  Sond 

Simons  ol  «L 

0.49  mm  Sond 

TIson 

004*179  mm  Sond 

TnSoki  of  ot 

1.03  mm  Sond 

TnSoM  of  ot 

120  mm  Sond 

” 

Tmkok)  of  at 

L40  mm  Sond 

TaoMii  of  a. 

2.20  mm  Sond 

• 

Oonoo 

10 


12 


14 


16 


18 


Fig.  5.  Comparison  of  predicted  and  observed  ranges  of  kd. 


Figure  B.4.3. 
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TABLE  II 

COMFAKISON  OF  COMPUTED  AND  MEASURED*  VALUES  OF  Ui  AND  2ao'** 


Run 

No. 

U 

ft/sec 

H 

m 

T*X10> 

ftVsec 

Ut 

(Equation  19) 
ft/min 

Ui 

(meas.) 

ft/min 

2«, 

(Equation  20) 
ft 

2<». 

(meas.) 

ft 

18 

2.1 

1.01 

0.04 

0.06 

0.06 

0.04 

28 

HQ 

3.7 

1.04 

musM 

0.03 

0.07 

0.12 

0.05 

IS 

119 

2.6 

1.05 

0.09 

0.14 

0.11 

0.09 

14 

1.60 

2.8 

0.58 

0.07 

0.07 

0.06 

0.03 

34 

1.64 

2.9 

0.54 

0.040 

0  08 

0.12 

0.06 

0.02 

16 

3,6 

0.168 

0.13 

0.13 

0.15 

0.16 

35 

1.80 

2.9 

0.53 

0.120 

0.19 

0.20 

0.08 

0.06 

17 

2.10 

3.5 

0.241 

0.19 

0.18 

0.16 

0.18 

33 

1.83 

3.3 

0.56 

0.161 

0.22 

0.18 

0.09 

0.08 

S 

2.21 

5.2 

0.93 

0.29 

0.17 

0.17 

0.20 

0.28 

10 

1.88 

3.9 

0.46 

0.23 

0.31 

0.24 

0.08 

0.12 

37 

2.54 

6.6 

1.11 

mwm 

0.25 

0.27 

0.28 

0.34 

36 

3.5 

0.55 

0.32 

0.37 

0.32 

0.10 

0.13 

6 

2.68 

5.8 

1.04 

0.44 

0.70 

0.45 

0.27 

0.31 

7 

2.14 

4.5 

0.59 

0.37 

0.35 

0.44 

0.12 

0.19 

38 

2.78 

6.3 

1.20 

0.52 

0.58 

0.27 

0.30 

11 

7.4 

0.92 

0.98 

0.46 

0.60 

0.25 

0.32 

8 

2.46 

5.9 

0.57 

1.03 

0.92 

0.59 

0.14 

0.17 

12 

3.12 

4.8 

0.89 

2.16 

1.16 

1.16 

0.22 

0.23 

*  (3)  pp.  68-9. 

^  »»2.64.  1.30  ft/sec,  ^ 1,  mean  sand  si2C"0.93  mm. 


Table  B.4.1 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  (metric)  units 
is  presented  on  page  iv. 
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Hfitation 


Most  of  the  notation  in  this  paper  is  the  same  as  in  Kennedy's  1963  paper  (see 
B.3).  New  notation  is  as  follows: 

So  =  equilibrium  amplitude  of  bed  fonns 

Tb  s  average  sediment  discharge  in  migration  of  the  bed  forms 

Uc  =  critical  velocity  for  initiation  of  sediment  motion 

^  =  ratio  of  bed-load  discharge  to  suspended-load  discharge 
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B.5.  Hayashi  (1970) 


Basis  of  theory 

1 .  Hayashi  modified  Kennedy's  (1963)  theory  by  adopting  a  slightiy  different 
transport  relation,  which  takes  account  of  the  effects  of  bed  slope  on  sediment  discharge: 

T(x,t)  =  m^l  +  a  [(|)x(x  -  5,  -d.t)]  (1) 

Principal  results  and  observations 

1 .  Modified  criteria  for  the  conditions  of  occurrence  of  different  bed  forms 
were  obtained  (Fig.  B.5.1). 

2 .  Agreement  of  theoretical  and  observed  conditions  for  occurrence  of  different 
bed  forms  is  quite  satisfactory  (Fig.  B.5.2). 

3 .  No  new  relations  for  heights  or  wavelengths  of  bed  forms  were  developed. 

4.  Hayashi  recommended  1.5<C<3.  C  =  2  appears  to  yield  quite  good 
results. 

5 .  This  paper  points  up  the  importance  of  local  bed  slope  along  the  bed  forms 
on  their  formation  and  behavior. 
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Figure  B.5.1 


FIG.  8. -COMPARISON  OF  FREDICTEI}  REGIONS  OF  OCCURRENCE  OF  SAND 
WAVES  FOR  CASE  C  •  2.0  WITH  EXPERIMENTAL  RESULTS  SUMMARIZED 
BY  KENNEDY  (2) 


Figure  B.5.2. 
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Notation 


Most  of  the  notation  in  this  paper  is  the  same  as  in  Kennedy’s  1963  paper  (see 
B.3).  New  notation  is  as  follows: 


C 

F 

Fi 

F2 

Fa 

Fm 


a/(2g5AJ2] 

Fronde  number 

maximum.  F  for  formation  of  dunes 


maximum.  F  for  formation  of  antidunes 
minimum.  F  for  function  of  antidunes;  F? 


F2forC  =  0;F^  = 


coth  kd 
kd 


T(x,t)  =  local  sediment  discharge 


Note:  (F^,  F^)  =  4icd  tinrk'd  ^  +  2±  V  (C+2)2  .  8  C  tanh^  kd  ] 


a  =  dimensionless  coefficient 
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B.6.  Engelund  (1970) 


Basis  of  theory 

1.  Vorticity  transport  equation.  A  stream  function  is  used  to  specify  the 
motion.  (Hence,  a  2-D  analysis). 

2 .  Continuity  equations  for  suspended  load  and  total  load. 

3 .  Transport  formulas  for  bed  load  and  suspended  load  discharges. 

4.  Assumption  that  eddy  viscosity  is  uniform  over  depth  and  there  is  a  "slip 
velocity"  at  the  bed  between  the  flow  and  the  bed. 

5.  Perturbation  stability  analysis.  A  sinusoidal  perturbation  is  introduced  in 
bed  level  and  water  surface. 

6.  Boundary  conditions:  (a)  vanishing  vertical  fluxes  at  bed  and  water  surface; 
(b)  no  shear  stress  at  water  surface;  (c)  appropriate  description  of  bed  shear 
stress  in  terms  of  slip  velocity;  and  (d)  relation  between  bed  shear  stress  and 
sediment  concentration  at  the  bed. 

Principal  results  and  observations 

1 .  When  bed  load  discharge  is  neglected,  bed  instability  is  predicted  for  a 
range  of  Froude  numbers  around  and  above  unity  (antidunes).  Upper  limit 
is  given  by 
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=  coth  (kD)kD 


(1) 


(same  as  predicted  by  Reynolds  1965),  and  lower  limit  is  given  by 

F2  =  tanh(kD)/kD  (2) 

(see  Fig.  B.6.1). 

2.  When  bed  load  discharge  is  included,  bed  instability  is  predicted  in  the 
lower  range  of  Froude  numbers  as  well  (dunes)  (see  Fig.  B.6.2). 

3 .  Complete  solution  is  seen  in  Figure  B.6.3. 

4 .  No  prediction  of  dune  height  is  given. 


32 


Eigures 


kD  =  (ZnlDD 

Fiouhe  5.  Stability  diagram  for  the  case  of  negligible  bed  load,  corresponding  to  VIU,  =  17. 

Figure  B.6.1. 


33 


0 


0-5 


0 


15 


kD 

Fiottke  6.  Asymptotic  stability  boundaries.  EIxperiments  by  Guy,  Simons  &  Richardson 
(1966):  O.  d  =  0-19-0-47  mm;  and  by  Kennedy  (1961);  x  ,  d  s  0-23-0'95  mm. 


Hgure  B.6.2. 


34 


Fiouke  10.  Stability  diagram  for  the  complete  solution.  The  parameters  are 
VIU,  =  21  and  V,l(wF)  =  1. 


B.6.3. 
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Notation 


d 

= 

grain  size 

D 

= 

flow  depth 

F 

= 

Froude  number 

k 

= 

wave  number 

Uf 

= 

friction  velocity 

V 

= 

mean  velocity 

w 

fall  velocity 
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B.7.  Gill  (1971) 


Sediment  continuity 


dz  ^  aqs  9v  _  Q 
3t  3v  3x 


Flow  continuity 


q  =  v(h  -  z) 


Sediment-transport  relation 


qs*Cqso^;^-  ij 


Friction  relation 


^  ^2g 


Momentum  equation 


3v  3h  gv2  _ 
V  — +  g  —  +  -5° - =  0 

(h-z) 
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1 .  Fot  the  case  of  zero  hricdon  (f  =  O,  Co  ->  «>),  the  bed-fonn  celerity  is 


’d(l-F2)r  V 


and  the  bed-form  height  (see  Fig.  B.7.1)  is  given  by 


2an  I  X 


2 .  The  corresponding  results  if  friction  is  included  are 


d(l-F2)(l-^) 

t 


[■•i'rg] 


d(l  -  F2)  I  - 
, _ _A _ L 

-  Ti  1 

2na  1  T  7  — 

2  ^  dv, 


3.  General  expression,  without  specification  of  a  sediment-transport  relation, 
for  vb  is 
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(10) 


V  3qs 

d(l  -  F2)  av 

The  relation  was  earlier  obtained  by  Reynolds  (1965)  and  Gradowcyzk 
(1968). 

4.  If  dynamic  effects  are  neglected  (F  =  0),  (7)  yields,  for  n  =  3, 


which  is  idendcal  (except  for  a)  to  Yalin's  (1964;  see  Section  C.2)  Eq.  13. 


5 .  Gill  presented  no  verificadon  of  his  results. 
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Notation 


C 

Co 

D 

f 

F 

g 

n 

q 

qs 

qso 

Ss 

t 

V 

Vb 


Some  notation  is  defined  in  the  figure. 

coefficient  in  (3) 

=  Chezy  coefficient 

=  particle  diameter 

=  Darcy- Weisbach  friction  factor 

=  Froude  number 

=  gravity  acceleration 

=  exponent  in  (3) 

=  water  discharge  per  unit  width 

=  sediment  discharge  per  unit  width 

=  specific  gravity 

=  time 

=  mean  velocity 

=  celerity  of  bed  forms 


-  1)  D3 


a 

Y 

t 

Xc 


1  2 

dune-shape  factor  (a  =  for  triangular,  a  =  —  for  sinusoidal) 

^  7C 

specific  weight  of  fluid 


bed  shear  stress 


critical  bed  shear  stress  for  initiation  of  motion 
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B.8.  Fredsoc  (1974). 


Basis  of  theory 

1 .  Same  as  that  of  Engelund  (1970) 

2.  Influence  of  gravity  on  bedload  transport  is  accounted  for  in  transport 
reladon: 


Ob  =  8(0  +  pi  -  0.47) 


(1) 


Principal  results  and  observations 

1 .  Stability  criteria  are  somewhat  different  hrom  those  obtained  by  Engelund. 
(see  Fig.  B.8.1). 

2.  By  investigating  growth  of  dunes  using  second-order-approximation, 
theory  explains  asymmetry  of  dune  shape  (steepening  of  the  downstream 
faces). 

3 .  No  relation  for  dune  height  is  derived. 
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Eism 


kD 

Figuke  2.  Stability  diagrams.  ViU„  =  20.  (o)/i  =  0.  (6)  ;»=  0- 1.  £.xpenments  by  Guy  et  al. 
( 1966):  X  ,  antidunes:  •.  standing  waves;  0>  dunes;  LL,  lower  limit  to  sediment  transport. 


Figure  B.8.1. 
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Notation 


d  =  particle  diameter 

D  =  (tepth 

F  =  Froude  number 

g  =  gravity  constant 

I  =  local  bed  slope 

k  =  wave  number 

qb  s  bed-load  discharge 
s  =  specific  gravity  of  sediment 

u  =  coefficient  (r  0.1) 


6 

<l>b 

X 


X 

g(s-l)d 

-  Ob _ 

[(s-l)gd3]i/2 

bed  shear  stress 
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B.9.  Richards  (1980) 


Basis  of  theory 


1 .  Equations  of  motion 


,  3u  ap  a  d 

U  —  +  W  —  =  -  —  +  —  tx2  +  —  Txz  +  g  sin  a 
ax  dz  ax  ax  “  az 


(1) 


u 


dx 


dp  a 
— + — 

dz  dx 


a 

+  -  '^zz » 

az 


(2) 


2.  Fluid  continuity 


au^^ 

ax  az 


(3) 


3 .  Transformed  coordinate  system 


X*  =  X,  z*  =  z  -  hoexp(ikx)F(z*) 


where 

sinh  k(D-z*) 
^  ~  sinh  kD 


(4) 

(5) 


4.  Qosure  relations  expressing  turbulent  stresses  in  terms  of  turbulent  kinetic 
energy  and  its  dissipation  rate  (equations  not  included  here). 


5 .  Sediment-transport  relation 
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qb  =  C(to)™ 


(6) 


Sediment-condnuity  reiadon 


9qb 

dx 


/I 

(1-n) — 
dt 


Perturbed  bed-load  discharge  taken  to  be 


(7) 


(8) 


Linearizadon  of  equadons  to  first  order  in  bed  forms  steepness  hok. 

Exandnadon  of  maxima  of  bed-shear-stress  variadon  along  x,  and  thus  also 
of  qb  (through  (6)),  to  find  wavelength  that  has  fastest  growth  rate. 

results  and  observadons 

Finds  two  maxima  in  bed  shear  stress.  Argues  that  one  corresponds  to 
ripples,  the  other  to  dunes. 

For  ripples 


for  1 .4  <  P  <  2.9.  L  is  independent  of  depth.  Bed  is  always  stable  for  P  > 
2.9. 


3 .  Stability  limits  and  dominant  wavelength  for  dunes  given  in  Fig.  B.9.2. 


4.  From  examination  of  roughness  length  Zq  for  flow  over  moving  beds, 
concludes 


203d  <  Xr  <  4,050  d 


(10) 


for  ripples. 

5.  Theoretical  prediction  of  dune  wavelength  (see,  e.g..  Fig.  B.9.2)  are 
consistent  with  results  of  flume  experiments  cited  by  Allen  (1970) 

Xd  =  2jtD  (11) 

and  Allen's  (1970)  empirical  relation 

X<i=  1.16  D1-55(X<,,  Din  m)  (12) 


(see  Fig.  B.9.3). 
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1  10  100 
Mein  depth,  D  (m) 


Fiocre  1 1 .  Comparison  of  observations  of  the  group  mean  wavelength  of  dunes  from  flume,  river 
and  marine  environments  (taken  from  Jsokaon  1976)  with  the  model’s  results:  O,  observations; 
•  ,  z,  s  0-t  cm;  A,  2,  »  0-25  cm;  +,  **  «  O-Scm;  ■,  z*  >  l-Ocm. 


Figure  B.9.3.  Equation  (5.6)  is  (11);  Equation  5.7  is  (12). 


49 


fioiatioa 


Some  notation  is  defined  in  the  figures. 


C 

d 

g 

k 

L 

m 

n 

Qb 

s 

% 

U* 

U,W 


Zo 


coefficient 
bed-particle  diameter 
gravity  acceleration 
IkIL 

X  =  wavelength 
exponent 

porosity  of  bed  material 
bed-load  discharge 

perturbed  (by  bed  waves)  bed-load  discharge 
bed  shear  velocity 

mean  velocity  in  x,z  direction 

curvilinear  coordinate  system  such  that  z*  =  const  are  (approximately) 

streamlines 

roughness  length 


P 

^,d 

^xy*  dC. 
to 


coefficient  defined  by  (8) 
wavelength  of  ripples,  dunes 
stress  (usual  subscript  notation) 
bed  shear  stress 
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B.IO.  Fredsoe  (1982) 


This  analysis  will  be  described  in  greater  detail  because  of  its  particular  relevance  to 


the  present  project. 


Sediment  continuity 


do  dh 

^  =  -(l-n)  — 
dx  at 


Integration  of  (1)  yields  for  bed  forms  migrating  without  change  of  form 


q  =  qo  +  a(l-n)h 


2.  If  all  material  is  transported  as  bed  load  (qo  =  0),  bed-form  celerity  is 


Qtop 

"(l-n)H 


3.  From  (2)  and  (3), 


H  qtop 


4.  Sediment-transpon  relation 
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where 


<Db  =  -p-^^  r.9*  =  — - 

V(s-l)gd3  P(s-l)d 


To  account  for  local  effects  of  bed  slope,  Fredsoe  adopts 


Then  (4)  becomes 


<i>nref 


Distribution  of  shear  stress  along  upstream  face  of  dune  (inferred  from 
experimental  results) 


ffi) 

e*  (.  if  W 

tq)  ^  2D  ■  dJ 


where  f  is  obtained  from  experimental  data  on  the  distribution  of 
boundary  shear  downstream  from  a  rearward  facing  step  on  a  rigid  bed. 


6 .  Adopt  Meyer-Peter  formula  to  get  form  of  (5): 


<D  =  8^l 


-  9q  -  0c 
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Then  (5),  (6),  (7)  and  (8)  yield  a  differential  equation  for  h 


where  0*  is  given  by  (9). 


(11) 


Principal  results  and  conclusions 


1.  Shape  of  bed  forms  is  obtained  by  integrating  (11)  numerically,  with 
upstream  boundary  condition  0  =  0c  at  h  =  0.  At  bed-form  crest,  0*  =  0jQp 

on  right  side  of  (1 1).  Downstream  face  of  bed  form  is  taken  to  be  at  angle 
of  repose  of  bed  material.  Bed  stress  is  obtained  from 


V  D’ 

-^  =  6  +  2.5  in^  (12) 

Uf 


Typical  shapes  of  bed  forms  are  shown  in  Fig.  B.10.2;  bed-form  steepness 
is  shown  in  Fig.  B.10.3. 

2.  Bed-form  height  is  calculated  by  introducing  a  disturbance  at  the  crest  and 
examining  conditions  under  which  it  is  stable.  Result  is 


Theoretical  dune-height  variation  with  8*  and  d  is  shown  in  Fig.  B.10.4. 


3 .  Wavelength  relation  is  given  by 

[16Hqb  +  (16H  +  S)qsl 
(qb  +  qs) 


(14) 


Friction  relation  is  obtained  by  applying  Carnot's  relation,  as  used  by 
Engelund,  and  the  theoretical  relation  for  h/H,  (4).  Result  is 


0  =  0'  +  0“  =  0' 


1_H  H 
2  D  L 


(15) 


4.  Combining  friction  and  wavelength  relations,  (14)  and  (15),  gives  bed- form 
steepness  as  shown  in  Fig.  B.10.5. 
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RG.  1,— Oaflnitlon  Skrtch  of  Mlflrotlns  Duno* 


Figure  B.10.1. 


RG.  S.— ExampiM  on  Calcuiatod  Duno  Profllot  and  Variation  in  Bod  Shoar  along  Thoao 
Profiloa:  (a)  6%  -  0.057.  HID  >  0.04.  0,  «  0.05;  (M  6%  -  OJO.  HjD  >  0.27.  8.  >  0.05 


Figure  B.10.2. 


f 


FK3.  6.— Variation  In  Langtti-HalvM  Ratio  witli  Bad  Shoar  Strata  lor  OWarant  Grain 
SIfao 


Figure  B.10.3. 
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FIQ.  8.— Variation  In  Ouna  Haight— Wator  Dapth  Ratio  with  Bad  Shaar  Straaa  for  ON- 
iarant  Orain  Slaaa 


Figure  B.10.4. 


FIG.  10 — Variation  in  langth-Haight  Ratio  with  Total  Bod  Shaar  Straaa  for  DNfarant 
Grain  Sizaa  (21) 


Figure  B.10.5. 
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Notation 


The  author's  full  notation  list  is  reproduced  here. 

a  =  migntioii  velocity  of  dunes; 

Of  »  Founer-coefficient  for  perturbation  of  bed; 

Cf  -  drag  coefficient; 

Cyo  =  maximum  drag  coefficient; 
d  «  mean  grain  diameter, 

D  =  water  depth; 

D '  =  mean  boimdary  layer  thickness  along  dune; 

F  =  Froude  number, 
g  =  acceleration  of  gravity; 
h  =  local  dune  height; 
h  perturbation  of  dune; 

H  =  maximum  dune  height; 

H,  =  height  of  upstream  step; 
k  -  bed  roughtM»s; 
k,  =  wavenumber  of  bed  perturbation; 

L  =  dune  length; 
n  =  porosity; 

=  bed  load  sediment  traiuport; 
q,  ^  suspended  sediment  transport; 

-  constant; 

q,  =  total  sediment  transport; 
q^  =>  sediment  transport  at  dune  top; 

4  =  change  in  sediment  transport,  due  to  perturbation  of  dune; 
i  relative  density  of  sand  grains; 

I  ==  time; 

(/*  *  near  bed  flow  velocity; 

U}  =  effective  shear  velocity; 

H>  3  fall  velocity  of  suspended  sediment; 

X  =  coordinate  in  flow  direction; 

7  =  consunt,  defined  in  Eq.  21; 

8  ^  phase  lag; 

8*  »  displacement  thickness; 

8o  =  upstream  boundary  layer  thickness; 
c  -  eddy  viscosity; 

q  =  water  surface  undulations; 

6'  ==  effective  dimensionless  shear  stress; 

6*  ==  dimensionless  shear  stress  due  to  form  drag; 

0^  =  critical  dimensionless  shear  stress; 

6*  =  local  dimensionless  shear  stress; 

6^,^  =  dimensionless  shear  stress  at  top; 

0tiu  ~  maximum  dimensionless  shear  stress; 

Ot;  =  local  dimensionless  shear  stress  including  effect  of  gravity; 
i  =  imaginary  unit; 

)!.  =  constant,  defined  in  Eq.  7; 
p  =  density  of  water; 

=  dimensionless  sediment  transport; 

(hb  =  dimensionless  bed  load  transport; 

=  dimensionless  suspended  sediment  transport;  and 
t*  s  bed  shear  stress. 
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B.  1 1 .  Haque  and  Mahmood  (198S). 


This  paper  is  the  culminadon  of  a  series  of  publicadons  (cited  in  the  paper)  by  these 
authors  on  the  subject  of  ripple  and  dune  geometry.  The  rirst  of  these  publications  was 
Haque's  (1970)  M.S.  thesis,  under  Mahmood,  submitted  to  Colorado  State  University. 


Basis  of  theory 


1 .  The  flow  is  treated  in  two  parts  (see  Fig.  B.l  1.1):  an  attached-flow  zone 
along  the  upstream  slopes  of  the  bed  forms;  and  a  downstream  wake-flow 
zone  in  which  the  free  streamline  is  calculated  in  the  basis  of  potendal-flow 
theory. 


2.  Admissible  bed-form  shapes  obtained  by  analyzing  the  potendal  flow  over 
an  infinite  train  of  bed  forms  (Fig.  B.l  1.2),  to  find  profile  meedng  two 
requirements: 

(a)  There  is  a  standing  eddy  downstream  from  each  bed  form. 

(b)  The  velocity  gradient  along  the  bed-form  profile  is  condnuous  at  the 
separadon  point 

3 .  Sediment  condnuity 


Jy 

dx  3t 


(1) 


4. 


Sediment-transport  relation 


5 .  Equations  and  boundary  conditions  of  potential  flow. 


6.  Case  of  finite  depth  and  rotational  inviscid  flow  is  handled  using  finite- 
element  techniques. 


1 .  Mathematical  relation  obtained  for  upstream  faces  of  bed  forms  is  shown  in 
Fig.  B.l  1.3,  along  with  bed-profile  data  from  ACOP  Canals  (in  Pakistan). 
Results  imply  homologous  shape  of  all  ripple  and  dune  profiles. 


2 .  Upstream  face  of  bed  forms  is  about  two-thirds  of  their  total  length. 


3 .  For  case  of  infinite  depth,  bed-form  steepness  is  given  by 

a  _  0.2463 
L  m 


(3) 


4.  Finite-element  analysis  yields  results  shown  in  Figs.  B.l  1.3,  B.l  1.4,  and 
B.l  1.5  for  steepness  and  relative  depth  of  bed  forms.  Optimization  is 
obtained  by  minimizing  differences  between  values  of  left-hand  and  right- 
hand  sizes  of  (2). 
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rne 

STmtAHLime 


FIG.  1.— SdiMiatie  Skateti  for  RIppIo  or  Dun* 


Figure  B.11.1. 


FLOW 


y 


a  \ 


/-* 


FIG.  3.— <«)  infintto  Soquonco  of  Attamatlng  Wedgts  «nd  Eddim;  and  (b)  DaHnL 
lion  Stateh  for  Oodtonn  Shapo  Functtona 


Figure  B.U.2. 
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i^VALU£  OF  EXPONENT  m  mREUATtVE  DEPTH  D/L 

RG.  9.— Effect  of  Row  Oeptti  on  Optimal  Bedfonn  Steepness 


Figure  B.11.3. 


RELATIVE  DEPTH  D/L 


RG.  10.— ReMlonsfiip  between  Optimal  Bedfonn  Steepness  and  Relative  Depth 
tor  Different  Values  of  m 

Figure  B.11.4. 
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OPTfUAL  SLOPC 


FLOW  DEPTH  RELATIVE  TO  BEOEORM  LENGTH  d/L 

FIG.  11.— Compartaon  to  Obsarved  Data 


Figure  B.11.5 


*  ^ 


Notation 


The  author's  full  notation  list  is  reproduced  here. 


a 

c 

D 

d 

E 

h 

k 

L 

L, 

m 

N 

t,t' 

u 

“c 


x,y 

x',y' 

3 


height  of  the  separation  point  from  the  stagnahon  point, 

celerity  of  bedforms, 

depth  of  flow  a'  the  crest, 

average  depth  of  flows, 

error  between  functions  V  and  N, 

bedform  height, 

dimensionless  length  of  the  upstream  face, 

bedform  length, 

eddy  length, 

a  dimensionless  index, 

the  bedform  shape  function, 

volumetric  sediment-transport  rate  per  unit  width, 

sediment-transport  rate  at  the  crest, 

time  in  the  fixed  and  the  moving  frame  of  reference, 

magnitude  of  local  velocity, 

velocity  tangential  to  bed  at  the  crest, 

undisturbed  velocity  at  infinity, 

(u/uj",  normalized  velocity  function, 

coor^nates  of  a  generic  point  in  the  fixed  or  moving  frame 
of  reference, 

coordinates  of  a  generic  point  in  the  moving  frame  of  refer¬ 
ence,  and; 

a  constant  roughly  equal  to  0.8. 
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C.  Empirical  Methods.  The  complexity  of  river-flow  phenomena,  and  the 
difficulties  encountered  in  treating  them  via  the  formalism  of  mathematical  fluid  mechanics, 
have  prompted  resort  to  purely  empirical  methods  of  correlating  the  variables  of  interest 
and  importance  to  river  engineers.  This  line  of  attack  was  begun  in  the  19th  century,  by 
proponents  of  "regime  theory"  who  sought  reliable  design  guidelines  for  the  irrigation 
canals  of  India,  which  transported  flows  with  significant  sediment  concentrations,  and 
often  presented  major  maintenance  problems  arising  from  bed  scour  or  deposition,  or 
unstable  channel  alignment  (bank  erosion).  It  was  only  natural  that  dimensional  analysis 
and  related  approaches  be  utilized  also  to  predict  and  quantitatively  describe  riverine  bed- 
forms. 


Eight  of  the  principal  bed-form  papers  of  this  type  are  summarized  in  this  section. 
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Cl.  Garde  and  Albertson  (1959) 


Basis  of  theory 


1 .  Dimensional  analysis  showing  that  dune  steepness  is  a  function  of  two 
dimensionless  quantides: 


h^ff—Ea— 

^  L(Ys  -  Yf)d 


(1) 


Principal  results  and  observadons 

1 .  Funcdonal  reladonship  for  (1)  is  obtained  by  curve  fitting  using  laboratory 
and  field  data  (Fig.  C.  1.1  for  ripples;  Fig.  C.1.2  for  dunes). 

2 .  Qitcria  for  different  flow  regimes  are  established  (Fig.  C.  1 .3). 

3 .  No  separate  relations  for  dune  height  and  dune  length  result  from  this 
analysis. 
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Figure  C.1.3. 
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Motaosa 


d  =  size  of  bed  material 

D  =  depth  of  flow 

Fr  =  Froude  number  (=  V/V^ ) 

h  =  average  height  of  ripples  or  dunes 

V  =  mean  velocity  of  flow 

Yj,  =  specific  weights  of  sediment ,  water 

X,  =  average  length  of  dunes 

Xq  =  average  shear  stress  on  the  bed 
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C.2.  YaUn(1964) 


Basis  of  theory 


1 .  Dune  geometry  as  shown  in  Fig.  C.2.1. 


2.  Assumption  that 


'CB='Cofi 


(1) 


and 


'CBsXcr 


(2) 


3 .  Dimensional  analysis  yields 


(3) 


Principal  results  and  observations 

1 .  Combining  ( 1 )  and  (2)  yields 


A 

D 


(4) 


If  the  variation  in  Xq/Xct  is  a  result  of  the  variation  of  flow  depth,  d,  only 
(i.e.,  S  =  const),  then 
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(5) 


tp  d 

2.  Functional  relationship  for  (4)  is  obtained  by  curve  fitting  using  laboratory 
and  field  data  (see  Fig.  C.2.2): 


3 .  For  rough  turbulent  fiow,  X  is  unimportant,  and  (3)  is  reduced  to 


Fot  hydraulically  smooth  flow,  d/D  is  unimportant,  and  (3)  is  reduced  to 


A 

p= const 


(8) 


4.  Functional  relationships  for  (7)  and  (8)  are  obtained  by  curve  fitting  using 
laboratory  and  field  data(see  Fig.  C.2.3): 


(X  >  ~  20) 

(9) 

const  -  KXX) 

(X  <  ~  20) 

(10) 

Hence, 

A  r  5d 

(X  >  -  20) 

(11) 
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AslOOOD.  (X<~20) 


(12) 


(1 1)  applies  to  dunes,  and  (12)  to  ripples 
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FIG.  1 


Figure  C.2.1. 
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Notaaon 


D 

= 

grain-size 

d 

= 

depth  of  flow 

F 

= 

Froude  number 

f() 

= 

function  of 

S 

= 

slope  of  uniform  flow 

V 

average  velocity  of  flow 

V, 

= 

Vxo/p  =  shear  velocity 

X 

= 

v^Dp/p  =  grain-size  Reynolds  number 

Y 

= 

specific  weight  of  fluid 

A 

s 

sandAvaveheight 

6 

= 

angle  of  repose 

A 

= 

sand-wavelength 

a 

viscosity 

P 

= 

density  of  fluid 

to 

shear  stress  on  the  surface  of  the  bed  =  ySd 

XB 

= 

shear  stress  at  point  B(see  Fig.  C2.1) 

4>() 

= 

function  of 

Subscript  (cr)  =  critical  value  signifying  beginning  of  panicle  movemenL 


The  bar  over  the  symbols  signifying  values  corresponding  to  flow  over  a  flat  bed  at  t 
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C.3.  Ranga  Raju  and  Soni  (1976). 


Basis  of  theory 


1.  Assumption  that  all  bed  forms  are  two-dimensional  and  of  identical, 
triangular  form. 

2.  Dimensions  (or  scale)  of  bed  forms  are  assumed  to  be  directly  related  to  rate 
of  bed-load  transport. 

3 .  Equation  for  bed-load  transport  per  unit  width 

qB=^)-U»-(l-X)-75  (1) 


4 .  Empirical  equation  for  migration  speed  of  bedform 


Uw 

VP^ 


=  0.021 


(2) 


5. 


Empirical  relation  between  bed-load  transport  and  grain  shear  stress 


Sb 

Ys 


^  Yf  ^ 

,Ys  -  Yf, 


1/2,  1  X  1/2  ^  . 


in  which 


Yi^bS 

(Yj  -  Yfld 


(3) 


(4) 
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1 .  Combinting  (1 ),  (2)  and  (3)  yields 


5F^F2  =  f2(Xo) 


in  which 


'Ys-Yf 


Functional  relationships  for  (5)  and  (6)  are  obtained  by  cur/e  fitting  using 
laboratory  and  field  data  (see  Figs.  C.3.1  and  C.3.2,  and  Table  C.3.1). 
Relation  obtained  from  Fig.  C.3.1  is 


^Fj  •F2  =  6.5x  103(to)8/3 


Rjj  is  obtained  from 


I  '  dl/6 

U  =  —  Rjj  ^  S  and  nj  =  24^  *>1  m.k.s.  units 


3 .  The  authors  find  that  more  than  75%  of  all  the  data  fall  within  a  ±  50%  band 
around  proposed  curve.  They  also  find  that  less  than  50%  of  the  data  fall 
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within  a  ±  50%  band  around  the  curve  proposed  by  Yalin  (1964)  (see 
Section  C.2).  Comparison  is  shown  in  Figs.  C.3.3  and  C.3.4. 
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Table  1.  Sources  of  data  and  Range  of  Variables 


no.  investigator 

ref. 

overall  range  of  parameters 

^mm 

/f»ft 

5x10* 

£/ft/s 

h  ft 

Lft 

1.  Guy  et  al. 

PI 

0.19  -0.93 

0.29(^1.070 

0.15  -  6.50 

0.70-3.43 

0.003-0.650 

0.29-17.7 

2.  Vanoni  and 
Hwang 

[141 

0.206-0.230 

0.197-0.870 

0.455-  2.90 

0.57-1.83 

0.037-0.057 

0.38-0.73 

3.  Vanoni  and 
Brooks 

[13] 

0.137 

0.227-0.305 

0.390-  2.80 

0.77-1.49 

0.042-0.061 

i 

0.34-0.45 

4.  Williams 

[18] 

1.35 

0.223-0.376 

1.33  -10.88 

1.50-2.63 

0.042-0.167 

1.30-8.90 

5.  Laursen 

[101  ! 

0.10 

0.180-0.797 

0.43  -  1.86 

1.07-3.36 

0.063-0.110 

0.45-0.56 

6.  Barton  and 

Lin 

[2] 

0.18 

0.293-1.190 

0.44  -  2.10 

0.74-3.60 

0.051-0.116 

0.43-0.74 

7.  Tsubaki 

[131 

1.26 

0.310-1.530 

1.61  -  1.73 

1.90-2.30 

0.070-0.270 

3.48-5.17 

8.  Bharat  Singh 

[1] 

0.60 

0.120-0.412 

1.00  -  7.00 

0.84-2.66 

0.009-0.066 

0.48-1.92 

9.  Martinec 
(Luanice 

River) 

[11] 

2.4 

0.470-4.62 

0.36  -  0.68 

1.16-3.22 

0.075-0.983 

1.13-11.9 

Table  C3.1 
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Notation 


The  author's  full  notation  list  is  reproduced  here. 


Notations 


Depth  of  flow 

(ijo  =  median  sediment  size,  for  which 
50%  of  the  material,  by  weight,  is 
finer _ 

£//Vg_«» _ 

Ulb-^d 

V  Qf 

Acceleration  due  to  gravity 
Average  height  of  undulations 
Average  length  of  undulations 
Manning’s  roughness  coefficient  for 
plane  bed 

Hydraulic  radius  of  the  bed 
Hydraulic  radius  of  the  bed  corres¬ 
ponding  to  the  grains 
Water  surface  slope 
Mean  velocity  of  flow 


Shear  velocity  =  'JgR^S 
Velocity  of  bed  wave 
Porosity  of  sand  mass 
Dimensionless  transport  parameter  = 

ifcTfe)' 

Bed  load  transport  in  weight  per  unit 
time 

Specific  weight  of  the  fluid 
Specific  weight  of  the  sediment 
Dimensionless  shear  stress 
corresponding  to  grains 
corresponding  to  grains  = 
yfRi,SKy,-'yf)d 
Critical  shear  stress 
Average  shear  stress  on  the  bed 
Kinematic  viscosity  of  the  fluid 
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C.4.  Yalin  and  Karahan  (1979). 


1 .  Dimensional  analysis  showing  that  dune  steepness  5  (=  A/A)  is  a  function 
of  1),  Z,  and  X: 


5  =  f(Ti.Z,X)  (1) 


where 


1 .  Analysis  of  data  showing  that  when  X  >  25  (dunes),  5  is  independent  of  X: 

5  =  <|)(ti.Z)  (5) 

2.  Analysis  of  data  showing  that  when  Z  >  r  100, 5  is  independent  of  Z: 

5  =  ({)(ti)  (6) 
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3 .  Laboratory  experiments  are  carried  out  to  determine  5's  dependence  on  Z  for 
Z  <  ~  100.  Results  are  plotted  in  Figure  2. 

4 .  Empirical  relation  is  established: 

g 

- —  =  ^exp(l-4)  (7) 

®inax 

5max  ~  (0.0127e*^)x  (8) 

in  which  ^  =  \/x  and  x  =  T|-1.  Fit  is  shown  in  Figure  3. 
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no.  1<— VMiatioa  of  Ouna  StupwiM  ««Mi  Rttothf  TiaMin*  Fotaa  tor  Laif*  ValvM 
M2 


Figure  C.4.1. 


no.  2.— Vrnaitofi  M  Oww  SlMparaa  mmIi  !!•>••>*•  Traciiv*  Fan*  tor  Variouo  Ranpoa 
M  Z 


Figure  C.4.2. 
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FIG.  3. — Relation  batwaan  Nortnslizad  Ouna  Staapnaaa  and  Tractive  Force  in  Loga¬ 
rithmic  Syatam  ol  Coordinataa 


Figure  C.4.3. 
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Notation 


cr  =  subscript  that  signifies  the  beginning  of  sediment  transport  ("critical 

stage"); 

D  =  representative  size  of  the  bed  material  (usually  selected  as  Dso); 

e  =  base  of  the  natural  logarithm; 

h  =  flow  depth; 

V*  =  shear  velocity; 

X  =  v*DA'  =  grain  size  Reyt.olds  number; 

X  =:  dimensionless  excess  of  the  tractive  force; 

X  =  value  of  x  corresponding  to  Smax; 

Z  =  h/D  =  dimensionless  flow  depth; 

A  =  dune  height; 

6  =  A/A  =  dune  steepness; 

5max  =  maximum  value  of  5  corresponding  to  given  (range  oO  Z; 

Tj  =  'iQ/iXo)a  =  relative  tractive  force; 

A  =  dune  length; 

V  =  kinematic  viscosity; 

^  =  x/x  =  normalized  value  of  x; 

p  =  fluid  density;  and 

to  =  shear  stress  acting  on  bed 
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C.5  Jaeggi  (1984). 


1 .  Analysis  of  laboratory  and  field  data. 


Upper  limit  of  alternate  bar-formadon 


ti  =  2.93  In  tib  =  3.13Z3^^ 


(see  Fig.  C.5.1). 


Minimum  slope  of  channel  necessary  for  alternate-bar  formation  is 


12.9;^ 

dm 


where  M  varies  from  0.34  for  uniform  bed  material  to  0.7  for  widely  graded 


material  (see  Fig.  C.5. 2). 


3 .  Scour  depth  due  to  alternate  bar  formation  is 


S  -0.76Aab  -  0.15 
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'n.jHz,'"* 


FIG.  3.^rlterlon  o(  Alternate  Bar  Formation  (7),  Compared  with  Teat  Reeulta 


Figure  C.5.1. 
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O'  10»  0>  10* 

Z,  'B/d 


FIG.  5.— Diagram  Giving  Minimum  Slope  for  Bar  Formation  for  Channel  Geometry 
and  Bed  Material  (11).  (Also  Indicated  are  Areas  which  Characterize  Some  Swiss 
Rivers  which  All  Fall  In  the  Corresponding  Region  (with  or  without  Bar  Formation)] 


Figure  C.5.2. 
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Notation 


B 

di 

dm  = 

M 

Pi 

s  = 

S 


channel  width; 

grain  size  within  the  bed  material  mixture; 

Jdidp;  mean  grain  size; 

parameter  which  is  function  of  sediment  grading; 
proportion  of  grain  size; 
reladve  density  of  bed  material; 
scour  depth; 


T1 

TIB 

e 

®cr 

0B 

Vm 


6/dcr,  ratio  of  Shields  factor  to  critical  Shields  factor, 

0B/0cr,  ratio  of  Shields  factor  characterizing  bar  forming  "horizontal 
movement"  to  cridcal  Shields  factor; 
hJ/(s  -  l)dm.  Shields  factor; 
cridcal  Shields  factor; 

BJ/(s  -  l)dm.  Shields  factor  related  to  channel  width 
mean  flow  velocity; 
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C.6.  Ikeda  (1984). 


Basis  of  theory 


1 .  Regime  relation  for  meander  wavelength 


X  =  otB?  =  aXr 


(1) 


2 .  Linear  stability  analysis  yields  bar  wave  length  of  maximum  instability 


X  =  a2^frc 


(2) 


3 .  Dimensional  analysis  yields 


pgDS  pU  D  B 

ViD’  (Ps-P)gd’  P  ’  D’  d 


(3) 


Hb  U  pgDS  pUD  _B  m 

D  (Ps-P)gd’  p  ’  D>  dj 


(4) 


Principal  results  and  observations 


1 .  Functional  relationships  are  obtained  by  curve  fitting  using  laboratory  data: 


-vs 


for  F  <  0.8 


(5) 


B  ^  ( d  1  D  P 


(6) 
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X  _  /B\0.55  /Dy  0.45 

(dj 


and 


Ha_ 

D 


=  0.0442  ( I 


-0.45  /B\1-9 


/BN 

UJ  w 


(7) 


(8) 


Figs.  C.6.1,  C.6.2,  and  C.6.3  present  verification  of  (1),  (6),  and  (8), 
respectively. 
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Notation 


a 

B 

Cf 

D 

d 

F 

F 

G 

g 

Hb 

s 

u 


numerical  constant; 
channel  width; 

resistance  coefficient  -  gDSAJ^; 
mean  depth  of  flow; 
sediment  particle  size; 
functional  relationship  in  Eq.  5; 
Froude  number  -  U/V^; 
functional  relationship  in  Eq.  6; 
gravitational  acceleration; 
bar  height; 
water  surface  slope; 
mean  fluid  velocity 


a  =  numerical  constant 

|3  =  numerical  constant 

ps,p  =  mass  density  of  sediment  and  fluid,  respectively 
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C.7.  van  Rijn  (1984). 


Basis  of  theory 

1.  Assumption  that  dimensions  and  migration  speed  of  bed  forms  are 
determined  by  rate  of  bedload  transport 

2.  Assumption  that  bed-load  transport  is  described  by:  (a)  a  dimensionless 
particle  parameter  D»;  and  (b)  a  transpon-stage  parameter  T 

(u*<r) 


3 .  Two  equations  for  rate  of  bedload  transport  per  unit  width: 


qb  =  SbUbCb 

(3) 

qb  =  (1  -  p)aAud 

(4) 

Principal  results  and  observations 

1 .  Combining  (3)  and  (4)  yields 

A 

Ch  Ub  Sb  D«io 

(5) 

d 

(1  -  p)a  “4*^50  d 
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It  is  assumed  that 


Cb 


8b  uii 
’D50’  Ud 


F(D*,T) 


(6) 


so  that 


A 

d 


.  d..t1 


(7) 


and 


.  D 


** 


(8) 


2.  Functional  relationships  of  (7)  and  (8)  are  obtained  by  curve  fitting  using 
laboratory  and  field  data  (see  Figs.  C.7. 1  and  C.7.2). 


j  =  0.1 1  (1  -  e  -  0-5T)(25  -  T) 


(9) 


and 


-  =  0.015 


^^j°\l-c-0-5T)(25-T) 


(10) 


Transport  stage  parameter  T  is  calculated  with  u*  given  by 
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where 


C  =  18  log 


The  author  notes  that  (9)  and  (10)  yield 


X  =  7.3  d 


which  is  in  agreement  with  Yalin's  result:  X  =  27tcL 


to- 


CO- 


- »  transport  stags  paramatar.T 

RG.  3.— Bed'Form  Steepness 

Figure  C7.2. 
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c 

Cb 

D 

D* 

d 

g 

P 

Rb  = 

s  = 

T 

u  = 

Ub 

Ud 

u*  = 


U* 


u*<r  = 


a  = 

A 

5b  = 

\ 

V  = 

P 

Ps 


Chdzy-coefficient  related  to  grains; 
bed-load  concentration; 
particle  diameter, 
particle  parameter; 
depth; 

acceleration  of  gravity  (LT*^); 
porosity; 

bed-load  transport  per  unit  width  (L^T'*); 

hydraulic  radius  related  to  bed  (L); 

specific  density  =  Ps/p; 

transport  stage  parameter 

mean  flow  velocity; 

velocity  of  bed-load  particles 

migration  velocity  of  bed  forms; 

bed-shear  velocity; 

bed- shear  velocity  related  to  grains; 

cridcal  bed-shear  velocity  for  inidation  of  modon; 

shape  factor  of  bedform; 
bed-form  height; 
thickness  of  bed-load  layer, 
bed-form  length; 

Idnemadc  viscosity; 
density  of  fluid; 
density  of  sediment; 
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C.8.  Menduni  and  Paris  (1986). 


Basis  of  theory 


1 .  Flow  over  each  dune  is  assumed  to  be  unaffected  by  upstream  dunes. 


2.  Momentum  equation  applied  to  control  volume  around  dune  (Fig.  C.8.1), 
which  yields 


(t-XB)L’  =  (pi-p2)H 


(1) 


3 .  Dimensional  analysis 

Principal  results  and  observarions 


1.  By  replacing  the  right-hand-side  of  (1)  with  drag  force  relation  obtained 
from  dimensional  analysis,  (1)  reads 

2 

(t-XB)L’=C*HpY  (2) 


By  rearranging  (2), 


_H 

L' 


(3) 
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(Note  that  there  are  several  typing  errors  in  the  authors’  equations.  These 
errors  have  been  corrected  hwein.) 


Total  shear  stress  is  calculated  as 

X  =  pgSh 


(4) 


Effective  shear  stress  is  calculated  as 

C  =  2.5  log 


(5) 


(6) 


2.  Laboratory  and  field  data  (see  Table  C.8.1)  are  used  to  determine  C*  as  a 
function  of  Re*: 


C*=^+K2l0(-K3^e*) 


Ki  =  10  (5.04-0.088B/h) 

K2  =  35 

K3=  1000 
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FLUME  DATA  FIELD  DATA 


discharge 

[m3/s] 

0.017  - 

0.357 

0.121  - 

1.751 

slope 

0.000075  - 

0.0017 

0.00037  - 

0.0079 

flow  depth 

[m] 

0.05  - 

0.47 

0.2  - 

2.5 

grain  size 

[mm] 

0.3  - 

1.46 

0.31  - 

1.44 

dune  height 

[m] 

0.003  - 

0.125 

0.015  - 

0.82 

dune  length 

[m] 

0.49  - 

3.25 

0.64  - 

19.22 

channel  width 

fm] 

0.08  - 

2.44 

35.69  - 

200.0 

Table  C.8. 1.  Experimental  data  ranges 


3.  Maximum  dune  steepness  is  calculated  by  setting  tb/x  =  0  (Fig.  C8.2).  One 
curve  is  for  "wide  channel"  (B  =  «>)  and  one  for  "narrow  channel"  (B  =  0). 

4.  A  comparison  of  measured  and  computed  dune  steepness  is  presented  (Fig. 
C.8.3). 

5.  The  authors  note  that  absolute  maximum  dune  steepness  is  0.07  (Fig. 
C.8.2). 
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04 


aunp 


0.1 


Fl9vire  5.  Maximum  theorecical  dun*  3t**Dness 

Rgure  C.8.2. 


10 


15 


duM  stMoncss 


Figure  C.8.3. 
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Notation 

B 

C*  = 
D 

g 

h 

H 

L 

P 

Re*  = 
S 

u*  = 


P 

X 

V 


channel  width 
form  drag  coefficient 

mean  grain  diameter 

acceleration  due  to  gravity 

flow  depth 

dune  height 

wavelength 

pressure 

dune  Reynolds  number  -  u^HA' 

slope  of  water  surface 
shear  velocity 

density  of  fluid 
total  bed  shear  stress 
effective  bed  shear  stress 
kinematic  viscosity 
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D.  Staristical  Models.  Following  World  War  II,  and  especially  during  the  1950's 
and  1960's,  wind  waves  (water  waves  generated  by  wind)  received  extensive  research  and 
engineering  attention.  Several  problems,  including  planning  for  the  wartime  amphibious 
landings  the  construction  of  many  off-shore  structures,  and  the  quest  for  better  harbor 
design,  prompted  this  wave  research;  review  of  these  is  beyond  the  scope  of  this  report. 
Suffice  it  here  to  note  that  the  differences  between  real  wind  waves,  which  include  fairly 
wide  spectra  of  wavelengths  and  periods,  and  monochromatic  waves,  which  had  been  the 
subject  of  most  of  the  classical  wave  studies,  turned  out  to  be  quite  important.  At  about  the 
same  time,  the  problems  of  describing,  managing,  and  separating  electromagnetic  waves 
comprising  many  frequencies  came  to  the  fore,  especially  in  radio  communication  and 
related  fields.  The  techniques,  notably  time-series  and  spectral  analysis,  used  in  describing 
random  and  other  multi-frequency  electromagnetic  waves  found  direct  application  in  the 
study  of  wind  waves.  (A  fascinating  and  comprehensive  description  of  these  developments 
is  presented  by  Kinsman  (1965)  in  his  classical  book  Wind  Waves.)  In  the  1960’s  the 
decidedly  non-monochromatic  character  of  riverine  ripples  and  dunes  gained  increasing 
recognition,  so  it  was  only  natural  to  apply  the  techniques  of  electromagnetic- wave  wind- 
wave  analysis  to  the  description,  and  to  some  extent  the  analysis,  of  these  bed  forms. 

Four  of  the  numerous  papers  of  this  type  are  reviewed  here,  and  others  are  cited. 
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D.l.  Nordin  and  Algert  (1966). 


This  paper  is  the  first  to  apply  autocovariance  and  spectral-density  analysis  to  bed 
forms,  and  as  such  was  a  major  contribution  to  the  subject.  The  study  reported  in  this 
paper  was  the  Masters  thesis  research  of  Algert  (1965)  carried  out  under  the  supervision  of 
Nordin.  The  work  was  continued  by  Nordin  as  the  subject  of  his  Ph.D.  thesis,  which  was 
also  reported  in  a  U.S.G.S.  Professional  Paper  (Nordin  1971). 

Basis  of  theory 

1.  The  established  mathematical  techniques  used  in  autocovariance  and 
spectral-density  analysis.  These  are  outlined  in  the  p^er. 

2 .  The  estimate  of  the  autocovariance  function  for  discrete  data  given  by 

N  -  X 

XX(i).X(i+X)  (1) 

3 .  The  estimate  of  the  spectral  density  function  for  discrete  data  given  by 

M  -  1 

G(f)-Co  +  2  ^  (2) 

i  =  1 

where  the  variance  Co  is  given  by 
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(3) 


Q)  =  ^  Z  X(i)2 
i  =  1 

Note:  It  might  be  useful  to  review  briefly  some  of  the  characteristics  of 
autocovariance  and  spectral  density  functions.  First,  it  should  be  noted  that 
one  can  speak  indifferently  of  the  autocoirelation  and  autocovariance 
functions  because  both  exhibit  precisely  the  same  properties;  the  former  is 
simply  the  autocovariance  divided  by  the  variance,  Co.  so  that  the  initial 
ordinate  is  unity.  Second,  the  spectral-density  function  represents,  in 
general  terms,  a  method  of  assigning  to  any  frequency  or  wave  number  a 
measure  of  its  contribution  to  the  "content"  of  the  process  or,  more 

specifically,  G(0  df  represents  the  contribution  to  the  variance  from  the 

^  ^  df 

frequency  range  f  ±  y. 

4.  Bed  profiles  (bed  elevations  measured  at  discrete,  regularly  spaced  points 
along  the  channels)  from  seven  flows:  3  in  a  4-foot-wide  lab  channel;  2  in 
an  8-foot-wide  channel;  and  2  in  the  Rio  Grande  conveyance  channel.  The 
data  for  these  flows  are  summarized  in  Table  D.1.1.  Typical  bed  profiles 
are  shown  in  Fig.  D.1.1,  and  the  autocovariance  and  spectral  density 
functions  for  Run  2  (see  Table  D.1.1)  are  shown  in  Fig.  D.1.2. 

5.  From  inspection  of  covariance  functions,  the  authors  concluded  that  the 
dune  process  can  be  represented  by  a  second-order  autoregressive  scheme, 
or  a  Markov  second-order  linear  model; 

X(i)  =  aX(i  -  l)-pX(i -2)  +  ei  (4) 
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Physical  reasoning  justifying  this  selection  is  presented  by  the  authors. 


Principal  results  and  observations 

1 .  Typical  observed  and  computed  spectral  density  functions  are  shown  in  Fig. 
D.1.3. 

2.  The  significant  wave  height  (average  amplitude  of  highest  one-third  of  the 
waves)  were  found  to  be  related  to  the  variance  by 

Hi/3  =  3cJ^^  (5) 

(For  ocean  wind  waves,  the  coefficient  is  4.) 

3.  H  1/3  was  found  to  be  a  linear  function  of  unit  discharge,  q,  for  all  three 
channels. 

4.  The  first  three  values  of  the  covariance  functions  relate  well  to  q,  suggesting 
that  the  signficant  statistics  to  generate  the  models  which  approximate  the 
process  are  functions  of  simple  flow  parameters. 

5.  The  temporal  autocorrelations  (at  a  point)  were  similar  to  the  spatial  ones. 
Longer  components  were  found  to  move  more  slowly. 

6.  Nordin's  (1971)  continuation  of  this  work  showed  that  the  spatial- 
frequency  spectra  can  be  normalized  as 
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G‘  =  ^  =  0.0001 


.8  J 


(6) 


(see  Fig.  D.1.4).  The  equilibrium  range  of  the  spectrum  is  seen  in  this 
figure  to  end  at  about 


(7) 


which  corresponds  to  a  wave  number,  k  =  27tf,  given  by 


F2  = 


0.19 

kd 


(8) 


This  wave  number  likely  represents  the  longest,  highest  bed  waves. 
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Figures  and  tables 


OISTAMCC  OOWNSmCAM  FflOM  MCAOSOX.tN  rCCT 


FIG.  2.— BED  PROFILES  FOR  DUNES  IN  LABORATORY  FLUME 


Figure  D.Ll 


W«VE  LENGTH  IN  EEET 


FIG.  4.— AUTOCOVARIANCF.  AND  SPECTRAL  DENSITY 
FUNCTIONS  FOR  RUN  NUMBER  2 

Figure  D.1.2. 
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FIG.  8. -OBSERVED  AND  COMPUTED  SPECTRAL  DEN¬ 
SITY  FUNCTIONS;  RIO  GRANDE  CONVEYANCE  CHANNEL 
NEAR  BERNARDO 


Figure  D.1.3. 
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G  r.) 


30 


Ficuiie  30.  -  Dimcnsioniesf  fp«cira  fi»r  ihe  proceai  i  ')■ 


Figure  D.1.4. 
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TABLE  1. -SUMMARY  OF  HYDRAULIC  DATA 


Run 

number 

Unit 

discharge, 

In  cfs  per  ft 

Depth,  In 
feet 

Water 

surface 

slope 

Number  of 
dunes  used 

Average 

dune 

length. 

In  feet 

(1) 

(2) 

<3) 

(4) 

(5) 

(G) 

1 

0.607 

0.400 

0.0037 

213 

1.0) 

2 

0.864 

0.485 

0.0038 

235 

2.26 

3 

1.05 

0.580 

0.0044 

131 

2.72 

4 

1.28 

0.670 

0.001 30 

20 

4.73 

5 

2.22 

1.05 

0.00131 

IG 

8.G3 

G 

O.IO 

2. GO 

0.00058 

13 

21.G 

7 

10.3 

4.15 

0.00058 

12 

25.4 

Table  D.1.1 


TABLE  2.-SUMMARYOF  STATISTICS 


Run 

number 

(1) 

C 

o 

(2) 

a 

(3) 

0 

(4) 

2 

(5) 

R 

(6) 

Significant 
dune  height. 

In  feet 

(7) 

Dune  length 
of  m.aximuni 
variance, 
in  feet 
(8) 

1 

0.00072G 

0.000264 

0.797 

0.081 

2 

0.00119 

1.205 

0.G17 

0.000326 

0.852 

0.104 

3 

0.00228 

1.309 

EMU 

0.000483 

0.893 

0.144 

3.34 

4 

0  01333 

1.480 

O.GIl 

0.003062 

0.856 

0.346 

8.4 

5 

0.03052 

1.170 

0.4G4 

0.003782 

0.936 

0.525 

13.4 

6 

0.419 

1.270 

0.518 

0.0918 

0.887 

1.94 

50 

7 

0.555 

1.095 

0.397 

0.184 

0.822 

2.24 

40 

Table  D.  1.2. 
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Notation 


The  author’s  full  notation  list  is  reproduced  here. 


A  »  dune  amplitude,  in  feet; 

Cq  =  variance  or  in  square  feet; 

Cx  =  autocovariance  functionrin  square  feet; 

E(f)  =  expected  value  of  f; 

f  3  wave  number  of  frequency,  in  cycles  per  unit  lag,  cycles  per  foot,  or 
cycles  per  second; 

G(f)  =  observ^  spectral  density  function  from  data,  in  cubic  feet; 
g(f)  =  spectral  density  function  model,  in  cubic  feet; 

Hi/3  =  average  of  one-third  highest  dunes,  in  feet; 

L  =  total  distance  interval  considered,  in  feet; 

I  =  distance  along  the  bed  profile,  in  feet; 

M  =  maximum  lag  used ; 

N  =  number  of  data  points  in  a  profile; 
r2  =  multiple  correlation  between  X(f) ,  X(1  -  1),X(£  -  2); 
r|(  =  correlation  coefficient,  C|(/Cg; 

T  =  time  period  for  a  dune  to  pass  a  given  point  in  minutes; 

X(l)  =  the  process  of  bed  elevation  X  as  a  function  of  distance  downstream, 
I,  in  feet; 

X(t)  =  the  process  of  bed  elevation  X  as  a  function  of  time,  t.  in  feet; 

Z^(f)  =  a  random  process  with  orthogonal  increments; 
a  =  autoregressive  coefficient; 

3  3  autoregressive  coefficient; 
r  »  dune  length,  in  feet; 

-  random  part  of  the  process,  X(t) ,  in  feet; 

9  3  angle  of  repose  of  bed  material,  in  degrees; 

X  ^  lag  interval,  in  feet;  and 
=  variance  of  Cj,  in  square  feet. 

V  «  mean  velocity 
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D.2.  Hino(1968). 


Basis  of  theory 

1 .  Observed  from  data  of  Nordin  and  Algert  (1969)  and  Ashida  and  Tanaka 
(1967)  that  log-log  plots  of  sand- wave  spectra  display  an  "equilibrium" 
range  wherein 

S(k)~k-3  (1) 

2.  Corresponding  power  law  (Kolmogorov)  for  turbulence  is 

E(k)~k-5/3  (2) 

3 .  Fourier  transformations  of  autocorrelation  functions  for  ti  and  Ti’  yield 

=  (4) 

4.  Sti'q'  has  dimensions  k'l.  Slope  t|'  of  sand  waves  cannot  exceed  angle  of 
repose  of  bed  material.  Therefore 

STi'T^k  =  (271)2  (x(0)k  -1  (5) 

5.  (4)  and  (5)  yield 

ST,T,{k)  =  a«|»k-3  (6) 
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Principal  results  and  ohservarions 


1.  Fig.  D.2.1.  shows  typical  results. 


2.  The  lower  limit  of  the  range  of  the  "-3  power  law"  is  given  by  (see  Fig. 
D.2.2) 


lco-0.15/h 

(7) 

L-7h 

(8) 

3 .  Ripples  (as  opposed  to  dunes)  tend  to  conform  to  a  "-2  power  law". 

4 .  C^antification  of  a(<|>)  from  experimental  data  led  to 


StjtiOc)  =  2.8  X  10-4k*3  (9) 


as  relation  for  spectral  density  function  in  the  equilibrium  range. 
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Figures 


k  (cyc/cm) 

FicrRE  1.  Wovp-nuraber  spectra  of  sand  waves  compared  with  the  ‘  —  3  power  law,  (10). 
(Curves  are  replotted  from  the  graphs  by  Nordin  &  AJgert  (1906).) 


Figure  D.2.1. 


20 


Notation 


E(k) 

h 

k 

L 

S'qt\(k)  = 
S'q*i^'(k)  = 


turbulence  spectral  density  function 
flow  depth 

wave  number  =  ML  (not  liML) 
wavelength 

bed-wave  profile  spectral  density  function 
bed-wave  slope  spectral  density  function 


a((j>)  =  same  function  of  <{> 

(j»  =  angle  of  repose  of  bed  sediment 

\|r(k)  =  wind-wave  spectral  density  function 
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D.3.  Annambhotla,  Sayre  and  Livesey  (1972). 


Basis  of  study 

This  study  was,  in  many  respects,  a  continuation  and  extension  of  the  work  of 
Algert  and  Nordin  (see  Section  D.l)  on  the  application  of  time-series  or  spectral  analysis  of 
bed  forms.  Its  principal  new  contributions  were: 

1.  Application  of  zero-crossing  analysis  to  bed  forms.  This  involves 
determination  of  a  best-fit  (according  to  whatever  measure  of  fit  is  adopted) 
straight  or  monotonically  curved  line  to  the  smoothed  bed  profile,  and 
analysis  of  the  distances  between  successive  crossings  of  the  best-fit  line  by 
the  actual  bed  profile;  see  Fig.  D.3.1. 

2.  Statistical  analysis  of  detailed  bed-form  data  from  a  moderatley  large  river 
(the  Missouri  River). 

3.  Filtering  of  the  raw  bed-profile  data  to  eliminate  the  very  low-frequency 
(long  wavelength)  components  due  to  meanders,  spur  dikes,  etc.  (see  Fig. 
D.3.2).  The  autocorrelation  and  spectral  analyses  were  based  on  the  filtered 
data.  Figure  D.3.3  shows  typical  autocorrelation  and  spectral-density 
functions  for  the  unfiltered  and  tiltered  data. 

4.  Frequency-distribution  analysis  of  the:  bed-elevation  deviations  from  the 
mean  (best-fit)  line;  wavelength;  and  amplitudes  of  the  bed  waves. 
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Principal  results  and  observations 


1 .  The  hydraulic  data  for  the  flows  that  produced  the  bed  profiles  analyzed  are 
given  in  Table  D.3.1. 

2.  The  results  obtained  from  the  spectral  analysis  are  summarized  in  Table 
D.3.2. 

3.  The  results  obtained  from  the  zero-crossing  analysis  are  summarized  in 
Table  D.3.3. 

4.  Weighted-average  bed- wave  characteristics  are  summarized  in  Table  D.3.4. 
These  were  obtained  by  weighting  the  wave  heights  and  lengths  by:  "the 
ratio  of  the  product  of  the  partial  discharge  assignable  to  a  sailing  line  and 
the  subreach  length  to  the  product  of  the  total  river  discharge  and  the 
combined  lengths  of  all  sailing  lines”  (pp.  506-7). 

5 .  An  overview  of  the  results  of  the  study  is  presented  in  Fig.  D.3.4. 

6.  The  bed-form  friction  factor  was  found  to  correlate  well  with  a  modified 
relative  roughness  (see  Fig.  D.3.5). 

7.  The  author’s  conclusions  concisely  state  their  principal  findings,  and  are 
reproduced  here  (p.  508); 

1 .  The  zero-crossing  distances  and  amplitudes  analysis  method  is  better 
suited  than  the  spectral  analysis  method  for  obtaining  characteristic  wave 
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heights  and  lengths  from  river  data.  This  is  mainly  because  the  former 
method  requires  only  stationarity  of  the  mean,  whereas  the  latter  requires 
both  mean  and  covariance  stationarity.  Bed  profile  data  from  rivers,  even 
after  filtering,  are  apt  to  be  nonstationary  in  the  covariance. 

2 .  Bed  elevations  were  found  to  be  distributed  approximately  according 
to  the  normal  probability  law.  The  approximation  was  found  to  be  better  for 
the  filtered  than  for  the  unfiltered  data. 

3 .  Wave  lengths,  amplitudes  and  heights,  determined  from  the  filtered 
data,  were  all  found  to  be  distributed  in  approximate  accordance  with  the 
exponential  probability  law. 

4.  The  bed  form  friction  factor,  /",  was  found  to  decrease 
approximately  exponentially  with  increasing  values  of  the  modified  relative 
roughness,  Rl(eH),  asymptotically  approaching  zero  as  the  contribution  of 
the  grain  roughness  to  the  total  resistance  becomes  predominant.  The 
parameter  was  found  to  vary  with  R/(eH)  in  a  manner  that  is  roughly 
consistent  with  the  von  Karman-Prandtl  logarithmic  law  for  rough 
boundaries.  Variations  in  the  size  distribution,  shape  and  geometrical 
arrangement  of  bed  forms,  which  are  not  adequately  accounted  for  in  the 
parameter,  RI(eH),  are  likely  contributors  to  the  scatter. 

5.  The  roughness  of  the  bed  was  observed  to  vary  inversely  with  water 
temperature.  However,  changes  in  Lq,  H,  e,  REleH,  and/'  all  tended  to 
lead,  sometimes  by  as  much  as  2  months  to  3  months,  rather  than  to  lag  the 
change  in  temperature.  The  reason  for  this  is  not  understood,  although  it  is 
suspected  that  unidentified  factors  such  as  seasonal  variations  in  sediment 
size  and  supply  may  play  a  role. 
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Figures 


FIG.  l.-DEFINTTION  SKETCH  OF  CHANNEL  BED  PROFILE 


Figure  D.3.1. 


•  KflllCt  III)  — 

aisTancc  •  (■■)  ^ 


(KftHtl  «  (l«»  — 


FIG.  ;i.-RAW  AND  FILTERED  BED  PROFILE  DATA  FOR:  (a)  RUN  R22:  AND  (A)  RUN 
RS8 


Figure  D.3.2. 
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A(«) 

0  10  to  SO  ao  SO  u 


FIG.  -t.-AUTOCORRELATION  AND  SPECTRAL  DENSITY  FUNCTIONS  FOR  UN¬ 
FILTERED  DATA  FROM  RUN  U22 


>  (•) 

0  10  to  so  «0  so  00 


FIG.  5.-AUTOCORRELATION  AND  SPECTRAL  DENSITY  FUNCTIONS  FOR 
FILTERED  DATA  FROM  RUN  R22 


Figure  D.3.3. 
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Tiai  IM  HOMTHt 


1  K;.  'J.-VA1«IAT10N  Of  WATKII  TEMPEKATHUK,  ItUlir.lINKSS  CIIAIIACTEIUSTICS, 
OLD  I'OKM  I'UICTIUN  KACTOH,  AND  ItIVEH  UlbCIIAKUE  WITH  IIME  OF  YEAK, 
lUUt),  IN  MISSOURI  lUVEK  AT  OMAiiA 


Figure  D.3.4. 
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FIG.  8.-VAUIATtON  OF  BED  FORM  FRICTION  FACTOR  WITH  MODIFIED  RELA 
TIVE  ROUGHNESS 


Figure  D.3.5 


TABLE  l.-HYnilAULIC  CONUITIONS  FOR  MISSOURI  lUVEU  HUNS 


130 


Table  D.3.1. 
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Table  D.3.2. 


TABLE  3.-WAVE  LENGTHS  AND  AMPLITUDES  FOR  REPRESENTATIVE  RUNS 
OBTAINED  BY  ZERO-CROSSING  AND  AMPLITUDES  ANALYSIS 


Wave  Length 

Positive 

Amplitude 

Negative 

Amplitude 

Wave  Height, 

/7  ■  a,  ♦  |o_l. 
In  (eet 
(meters) 

Run 

number 

m 

In  feet 
(meters) 

■ 

In  (eet 
(meters) 

c,. 

Roughness 
concentration 
f  ■=  11/7^ 

(1) 

HI 

H 

(4) 

(6) 

(7) 

(8) 

(9) 

RIT 

36.65 

1.690 

0.287 

1.280 

0.251 

1.574 

0.538 

0.0147 

(11.20) 

(0.088) 

(0.076) 

(0.164) 

R18 

43.99 

0.929 

0.473 

0.765 

0.582 

1.325 

1.015 

0.0231 

(13.40) 

(0.144) 

(0.177) 

(0.310) 

R19 

43.66 

0.739 

0.662 

0.741 

0.525 

0.873 

1.130 

0.0258 

(13.32) 

(0.202) 

(0.160) 

(0.345) 

R20 

21.60 

1.335 

0.300 

1.031 

0.305 

1.578 

0.0280 

(6.60) 

(0.093) 

(0.185) 

R21 

44.67 

0.794 

0.396 

1.131 

0.344 

1.373 

0.740 

0.0165 

(13.62) 

(0.121) 

(0.105) 

(0.226) 

R22 

32.57 

0.800 

0.557 

0.870 

0.418 

0.953 

0.976 

0.0300 

(9.93) 

(0.170) 

(0.128) 

(0.298) 

R23 

60.42 

0.712 

0.722 

0.692 

0.589 

2.150 

1.188 

0.0196 

(18.40) 

(0.220) 

(0.180) 

(0.362) 

Table  D.3.3. 


TABLE  4,-AVEaAGED  PROPERTIES  OF  BED  FORMS 


— 

Welidtlwl  Average  Wave  Charaeit-rlstlcs 

Cbic 

Description  of  rlvvr  bed 

L  , 

«> 

Iff 

In  feet 

In  feet 

.  -•  ft/L 

R/tU 

(1) 

(meiersi 

(metersl 

0 

(31 

(3) 

(41 

(5) 

(6) 

(7) 

Kf'US 

.1/25/66 

Dunes,  liars,  plane 

IS.O 

0.353 

0.0169 

3,910 

i<ii>  in 

1/8/66 

<I.S8) 

(0.077) 

Oars 

30.8 

0.3G2 

0.012G 

3,130 

4/20/68 

<€.36) 

(0.080) 

Uincs 

33.0 

0.845 

0.0353 

812 

lltl'IIIO 

5/13/C8 

(10.26) 

(0.258) 

IXines,  bars 

11.2 

0.827 

0.0187 

894 

HiT-lllU 

tl3.S0) 

(0.252) 

Dunes,  bars 

41.3 

0.872 

0.0311 

772 

iiju-  iij;i 

6/20/68 

(12.01) 

(0.266) 

Dunes,  bars 

39.3 

0.921 

0.0335 

045 

Mjii-  iUO 

r/i5/ii8 

(12.00) 

(0.383) 

[Xincs,  bars 

<2.5 

0.C36 

0.0150 

1,310 

iiai-uj7 

8/8/68 

(12.981 

(0.191) 

Transition,  dunes. 

36.0 

0.115 

0.0115 

3,330 

8/28/68 

bars 

(11.001 

(0.127) 

Transition,  plane 

31.2 

0.175 

0.0U512 

12,800 

2/27/68 

(10.  ID 

(0.053) 

KiU*  UlS 

Transition,  plane, 

(hincs 

35.1 

17.651 

0.163 

(0.050) 

0.00650 

10,G00 

UMt-lll'J 

10/10/08 

Plane,  transition 

10.1 

o.in 

0.0113 

4.950 

Krt0-K52 

10/24/08 

<3.081 

(0.011) 

Plane,  transition. 

21.3 

0.092 

0.00432 

3G.OOO 

ur,j-uo5 

long  liars 

(6.501 

(0.028) 

11/21/08 

Plane 

8.3 

0.089 

0.0108 

8,830 

(2.53) 

(0.037) 

Table  D.3.4. 
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Notation 


The  authoi's  full  notation  list  is  reproduced  here. 


c,, 

e 

f 

/' 

/" 

.K 

n 

}V 

H(f) 

!//(/)! 

h 

L 

1-0 

M 

Mom,. 

A' 

"o 

RyM 

S 

S' 

S" 

If 

X 

r 

r(/) 

r*o) 

v(a) 

v'(v) 

|v(v)| 

a 

k 

I' 

U  2.7/ 


positive  amplitude; 

negative  amplitude; 

coefficient  of  variation; 

roughness  concentration,  defined  as  H/Lg‘, 

Darcy-Weisbach  friction  factor;  also  frequency  or  wave  number; 
component  of  friction  factor  associated  with  sand  grain  roughness; 
component  of  friction  factor  associated  with  form  roughness; 
physically  realizable  one-sided  spectral  density  function; 
gravitational  acceleration; 

dune  height  defined  as  vertical  distance  between  maximum  and 

minimum  elevation  between  successive  zero-crossings; 

index  of  dune  height  (equal  to  2<Ty); 

frequency  response  function; 

gain  of  filter; 

arbitrary  wave  height; 

equivalent  sand  roughness; 

wave  length; 

wave  length  obtained  from  zero-crossings  analysis; 
maximum  number  of  lags  in  spectral  estimation; 

1th  moment  of  spectral  density  function; 

number  of  data  in  sample  of  bed  profile; 

number  of  zero-crossings  in  given  length,  X,  of  sample; 

autocovariance  function  for  y(.v); 

energy  slope; 

component  of  energy  slope  associated  with  sand  grain  roughness; 
component  of  energy  slope  associated  with  form  roughness; 
mean  flow  velocity; 

length  of  sample  record;  also  a  random  variable; 
distmee  along  length  of  channel  or  sample  record; 

Fourier  transform  of  y(.\ ); 

Fourier  transform  of  ,v'(x); 
bed  elevation  at  distance  x; 

bed  elevation  at  distance  x  for  raw  data  before  filtering; 

stochastic  process  defining  bed  elevation; 

standardized  form  of  random  variable  yCv); 

parameter  of  filter; 

spectral  width; 

sampling  interval; 

von  Karman’s  constant; 

lag  distance; 

kinematic  viscosity; 

standard  deviation  of  v;  and 

circular  frequency. 
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D.4.  Jain  and  Kennedy  (1974). 


1 .  Bed  profile  is  expressed  in  terms  of  its  Fourier  components. 


T](x,t)  =  i  B(k,t)  e**  dx 


Velocity  potential  for  flow  shown  in  Fig.  D.4.1  over  wavy  bed  given  by  (1) 


<|>(x,y,t)  =  iU  J  JB(k,t)e‘*“  dx 


where 


V  cosh  ky  -f  F^kd  sinh  ky 
F^kd  cosh  kd  -  sinh  kd 


3 .  Sediment  continuity 


t1i  +  Tx  =  0 


Note:  subscripts  denote  partial  diffemtiation. 


4 .  Transport  relation  of  Hayashi  ( 1970) 


T(x,t)  =  m(l  +  aTix(x,t)l[(U  -  Uc)  +  <Dx(x  - 13.  -  d.t)]" 
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5. 


(4)  and  (5)  yield  T|t(x,b),  which  is  Fourier  transfomied  to  yield 


Bt(k,t)  -  Tk2[a  -  iniJiexp(-  ikp)]B(k.t)  =  0  (6) 

after  substitution  of  (5)  and  (6). 

6.  Solution  of  (6)  is 

B(k,t)  =  B(k,0)exp[Tk2{a  -  iniJiexp(-  ikP))]  (7) 

7.  Spectrum  of  bed  profile  is  obtained  by  ensemble  averaging  the  product  of 
(7)  and  its  complex  conjugate 

<D(k,t)  =  <I>(k,0)exp[2Tk2(a  -  niJisinkp)t]  (8) 

8.  Speed  c  of  contribution  to  <1>  with  wave  number  k  is  found  from  Eqs.  1  and 
7  to  be 

c  =  TnikJ  icoskP  (9) 


9 .  Normalized  growth  rate  of  spectrum  is 


Ot(k.O)  d2 
<I>(k,0)  2Ta 


=  k2d2^1  -  ’^JisinkpjsF 


(10) 


which  has  a  singularity  at  the  Airy  celerity 
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=  tanh  kd/kd 


(11) 


Component  with  wave  number  given  by  (6)  develops  linearly  with  time 
even  on  a  flat  bed: 


■  =,  U  Ak2  w  ox 
^  (U  -  Uc)  sinh  kd 


10.  Normalized  growth  rate  given  by  (10)  is  shown  in  Fig.  D.4.2. 


11.  A  "conservation  of  variance"  principle  is  formulated  as  (see  Fig.  D.4.3) 


<I)(k,i)dk  =  A(k,t)dk  +  —  i  f  <I>(k,t)kc(k,t)dk  -  f  0(k,t)kc(k,t)dk 

2^1  Ik  Lk 


—  I  /  <I>(k,t)kc(k,t)dk  Jdk 


1 2.  For  sufficiently  large  kd,  Ji  approaches  infinity  and  Eq.  (9)  becomes 


c-Tnik 


1 3 .  For  steady-state  conditions  (equilibrium  range),  ( 1 3)  becomes 


^  k  J  Ocdk  =  0 
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which,  with  (14),  yields 


0~lr3 


(16) 


Principal  results  and  observations 

1 .  Spectra  of  bed  forms  developing  on  an  initially  flat  bed  exhibit  two  peaks, 
as  predicted  by  (10)  (see  Fig.  D.4.3).  One  peak  is  at  the  Airy  celerity  given 
by  (11). 

2 .  Authoi's  explanation  of  evaluation  of  spectra  is  as  follows: 

At  small  times  the  spectra  of  bed  forms  developing  on  an  initially  flat  bed 
are  characterized  by  two  peaks.  One  peak,  generally  that  at  the  lower  spatial 
frequency,  traces  its  origin  to  the  velocity-field  penurbation  and 
accompanying  pattern  of  differential  deposition  and  scour  on  the  bed 
produced  by  a  small  amplitude,  stationary  surface  wave;  i.e.  a  wave  train 
moving  relative  to  the  fluid  with  velocity  just  equal  in  magnitude  but 
opposite  in  direction  to  the  mean  flow  velocity.  It  is  this  equality  that 
determines  the  frequency  of  one  spectral  peak.  The  second  peak 
corresponds  to  bed  waves  resulting  from  the  inherent  instability  of  an 
interface  between  an  erodible  bed  and  a  turbulent  flow.  Over  wide  ranges 
of  flow  conditions  any  small  initial  disturbance  on  an  otherwise  flat  bed  will 
produce  a  perturbation  of  the  velocity  field  and  hence  also  of  the  sediment- 
transpon  distribution,  giving  rise  to  a  spatial  pattern  of  scour  and  deposition 
that  produces  bed  waves.  This  instability  mechanism  is  not  dependent  on 
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the  proximity  of  a  free  surface;  it  also  produces,  for  example,  ripples  in  very 
deep  tidal  flows  and  aeolian  ripples  on  wind-swept  sand  deposits.  The 
spatial  frequency  at  which  it  is  centered  is  a  function  of  the  flow  and 
transport  characteristics,  reflected  in  ni/a,  Ji  and  |Vd,  and  can  be  expected 
to  occur  at  the  value  of  kd  where  T,  given  by  (10)  and  illustrated  in  figure  2, 
has  a  continuous  maximum.  Equation  (9)  indicates  that  the  two  different 
families  of  waves  will  move  with  different  speeds,  and  this  frequency 
dispersion  probably  accounts  for  the  shift  of  the  normalized  spectra — 
towards  the  lower  frequencies  as  the  bed  configurations  develop  and 
mature. 

This  "variance-cascade"  process,  involving  interaction  between  the  water- 
surface  (Airy)  wave  generated  bed  forms  and  the  bed  forms  produced  just 
by  the  bed-flow  interaction  was  verified  by  the  experiments  of  Nakagawa 
and  Tsujimoto  (1984). 

3 .  Normalized  (by  do^F^/^)  spectra  conform  to  (16)  (see  Fig.  D.4.4). 

4.  Nonnalization  by  suggests 

L2/d  s  (M0/M2)  ~  F5/3  (17) 

which  is  supported  by  Jain's  data  (see  Fig.  D.4.5).  The  data  given  in  Fig. 
D.4.5  yield 

F5/3=^^  =  ^  (18) 

2jt  d  k2d  ^  ^ 
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Nordin's  data  confonn  better  to 


Ud-F2  (19) 

The  field  data  of  Annambhotla  et  al.  (1972)  do  not  conform  to  the  relation 
suggested  by  Fig.  D.4.5. 

5 .  Examination  of  bed-wave  spectra  of  several  authors  suggest  to  be 
preferred  in  the  spectrum  normalization  factor. 
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Ficnir  j.  Schomatic  rcprcscniation  of  differential  n-liuioii  for 
bed- wove  vorianec.  equation  ( 1!’)- 


Figure  D.4.3. 


10- 


Figuuk  0.  Non-diinensionol  plot  of  Jain’s  (1971)  spectra  data  illustrating  occurrence  o 
minu8•tllree-po^vcr  law  at  higher  frequencies. 

Uirtjs)  I  410  1  C’3  1  145  1  •J(i2  1-330  1  442  1  350  1  500  I  310  I  320  1  HO 

d  (ft)  0-329  0-331  0-248  0-253  0-267  0-240  0-437  0-402  0  335  0-353  0-417 


Figure  D.4.4. 


Figure  T.  Rulation  between  F  and  norrnalizcd  wavclcnsth  for  .fain  s  (11)71) 
and  Nordin's  (1971)  data.  #.  Jam's  data;  C-  Xordin's  data. 


Figure  D.4.5. 
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Notation 


Some  of  the  notation  is  defined  in  Fig.  D.4. 1. 


B 

c(k) 

F 

g 

i  = 

J 
Jl 
J'l 

k 

L 

m,n  = 
M0M2  = 
ni 

t  = 

T 

t 

Uc 


function  defined  by  (1) 

celerity  of  bed-wave  component  with  wave  number  k 
Froude  number  =  U/Vgd^ 
gravity  constant 

function  defined  by  (3) 
function  defined  in  (6) 

-J(k'.-d) 

wave  number  =  2k/L 
wavelength 

coefficient,  exponent  in  sediment  transport  relation,  (5) 
zeroth,  second  moment  of  spectral  density  function 
nU/(U-Uc) 
time 

local  volumetric  sediment  discharge  per  unit  width 
average  volumetric  sediment  discharge  per  unit  width 
critical  velocity  for  initiation  of  motion 


a  =  coefficient  of  bed-slope  term  in  (5) 

|3  =  sediment-transport  lag  distance 

r  =  normalized  spectral  growth  rate;  (10) 

<j)  =  velocity  potential 

(1>  =  spectral  density  function 
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III.  CURRENT  STATUS  AND  PROSPECTS 


A.  Introduction.  Section  B  of  this  chapter  will  pr^is  what  presently  is  known 
about  the  formation,  behavior  and  characteristics  of  alluvial  bed  forms,  and  the  principal 
deficiencies  in  our  knowledge  about  them.  Section  C  is  given  over  to  a  (somewhat  philo¬ 
sophical)  discussion  of  the  stumbling  blocks  that  arise  in  formulating  theories  of  bed-form 
behavior.  The  focus  of  this  chapter  is  very  long  riverine  dunes  (so-cailed  megadunes). 

B.  What  Is  Known.  The  following  features  of  riverine  bed  forms  and  their 
behavior  may  be  considered  as  fairly  well  to  finnly  established: 

1 .  Bed  forms  are  the  result  of  an  instability  that  is  inherent,  under  most  flow 
and  sediment  conditions,  in  the  interaction  between  the  primary  shear  flow 
and  the  bed,  or  between  velocities  induced  by  standing  (stationary)  surface 
waves  and  the  bed  (Exner  1925;  Anderson  1953;  Kennedy  1963,  1964, 
1969;  Reynolds  1965,  Hayashi  1970;  Engelund  1970;  Fredsoe  1974, 
Richards  1980,  and  others).  The  underlying  mechanisms  of  the  instability 
are  illustrated  in  Fig.  B.4.1.  Note  that  a  free  surface  is  not  required  for  the 
formation  of  bed  forms  (except  antidunes).  The  principal  processes 
involved  in  this  instability  are  as  follow: 

a.  Something  gives  rise  to  a  disturbance  (small  bump)  on  an  otherwise 
initially  flat  sediment  bed;  or,  in  the  case  of  free-surface  flows, 
produces  a  standing  water-surface  wave. 

b.  The  bed  disturbance  perturbs  the  flow  field  near  the  bed,  and 
thereby  perturbs  the  local  sediment  transport  rate. 
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c.  If  the  distribtuion  of  local  sediment-transport  rate  is  perturbed  so  as 
to  produce  a  pattern  of  scour  and  depositon  that  increases  the 
amplitude  of  the  bed  disturbance,  ripples,  dunes,  or  antidunes 
result,  depending  on  the  flow  and  sediment  characteristics. 

d.  If  the  distribution  of  local  sediment-transport  rate  is  perturbed  so  as 
to  produce  a  pattern  of  scour  and  deposition  that  diminishes  the 
amplitude  of  the  inidal  bed  disturbance,  a  flat  bed  results. 

e.  The  bed  forms  migrate  downstream  (ripples  and  dunes)  if  scour 
occurs  on  their  upstream  (stoss)  slopes  and  deposition  on  their 
downstream  (lea)  slopes.  Upstream  migration  (antidunes)  results  in 
the  opposite  case. 

2.  Over  a  fairly  wide  range  of  flow  and  sediment  conditions,  intermediate  to 
those  that  produce  ripples/dunes  and  antidunes,  the  flow-bed  interface  is 
stable,  and  remains  flat.  Any  induced  bed  disturbance  diminishes  in 
amplitude  and  disappears. 

3 .  The  bed  configuration  has  a  dramatic  effect  on  riverine  friction  factors.  The 
Darcy-Weisbach  friction  factor  often  changes  by  a  factor  of  ten  or  greater 
between  flat-bed  and  ripple/dune  conditions  (Vanoni  1975,  pp.  1 14-118). 

4.  There  are  no  generally  accepted  definitions  or  criteria  for  distinguishing 
between  ripples  and  dunes.  Their  behaviors  are  quite  similar.  Dunes  are 
longer  and  higher  than  ripples,  typically  by  a  factor  of  five  and  greater. 
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Ripples  often  occur  in  the  stoss  slopes  of  dunes.  There  does  not  appear  to 
be  a  continuous  spectrum  of  wavelengths  between  ripples  and  dunes;  on 
plots  involving  bed-form  wavelength,  dunes  appear  as  an  isolated  group  of 
points  (e.g.,  the  points  in  the  lower-left  of  Fig.  B.4.3). 

5.  There  is  no  readily  discemable  difference  between  the  spectral  density 
functions  of  ripples  and  dunes,  except  perhaps  for  their  equilibrium  -  range 
slopes  (see  item  13,  below). 

6.  A  phase  shift  between  the  local  bed- wave  displacement  and  the  local 
sediment  discharge  is  essendal  to  the  initiation  and  growth  of  bed  forms  (see 
references  cited  in  item  1.  above).  The  factors  that  have  been  identified  as 
contributing  to  this  shift,  and  their  relative  importances,  are  (Kennedy  1978) 
(note:  1  =  very  important  ->  4  =  not  important): 

a.  Phase  shift  between  bed  shear  stress  and  bed  waves  (i.e.,  bed 
displacement  and/or  slope),  due  to  effects  of  nonuniform  streamwise 
pressure  gradient  on  the  velocity  distribution.  (1  for  ripples,  2  for 
dunes,  3  for  antidunes). 

b.  Phase  shift  between  bed  waves  and  suspended- load  transport  rate, 
caused  by  time  required  for  sediment  to  be  entrained  into  an 
accelerating  flow  and  to  settle  out  of  a  decelerating  flow.  (1  for 
antidunes,  2  for  dunes,  3  for  ripples). 
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c.  Gravitational  force,  which  impedes  particle  motion  up  stoss  slopes 
and  aids  it  down  lee  slopes.  (1  for  ripples,  2  for  dunes,  3  for 
antidunes). 

d.  Flow  separation  from  downstream  slopes.  (1  for  fully  developed 
ripples  and  dunes,  4  for  antidunes). 

e.  Particle  inertia.  (3  for  antidunes;  4  for  others). 

f.  Phase  shift  between  local  flow  shear  and  local  turbulence 
characteristics;  i.e.,  the  "history"  effect  on  turbulence.  (2  for 
antidunes,  3  for  dunes,  4  for  ripples). 

g.  Phase  shift  between  bed  and  water-surface  waves,  and  hence  also 
between  bed  waves  and  local  mean  depth  and  velocity.  (2  for 
antidunes,  3  for  dunes,  unimporant  for  ripples). 

7 .  Practically  all  of  the  successful  analytical  models  developed  to  date  for  bed 

forms  are  based  on  small-perturbation  theory.  This  involves: 

a.  Deformation  (mathematically)  of  the  sediment  bed  into  a  small- 
amplitude,  slowly  moving  wave  (usually  a  sinusoid). 

b .  Calculation  of  the  flow  field  over  the  deformed  bed. 
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c .  Calculation  of  the  local  sediment  discharge,  along  the  deformed  bed, 
produced  by  the  perturbed  flow  field. 

d.  Application  of  the  sediment-continuity  equation  (Exner  relation)  to 
calculate  the  subsequent  defonnation  of  the  bed. 

e.  Interpretation  of  the  calculated  bed  profile  to  explain  the  occurrence 
and  behavior  of  different  types  of  bed  fonns. 

8.  The  analyses  outlined  in  item  3,  above,  are  limited  to  small-amplitude  bed 
forms,  and  cannot  explain  several  important  features  of  fully  developed  bed 
forms.  In  particular,  they  cannot  predict  the  behavior  and  properties  of 
mature,  fully  developed  bed  forms  which  are  dominated  by  nonlinear 
effects. 

9 .  Three  principal  items  are  forthccnning  from  most  of  these  small-perturbation 
analyses: 

a.  The  dominant  wavelength,  calculated  as  the  wavelength  with  the 
fastest  amplitude  growth  rate. 

b.  Bed- wave  celerity  as  a  function  of  wavelength,  the  flow 
characteristics,  and  the  sediment  and  sediment-transport  propenies 
of  the  flow. 

c.  Conditions  for  occurrence  of  the  different  types  of  bed  forms  (see, 
e.g..  Fig.  B.4.2).  These  predictions  invariably  involve  parameters 
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involved  with  the  bed-elevation-transpon-rate  phase  shift,  and  are 
difHcult  to  apply  to  predict  occurrence  of  different  types  of  bed 
forms. 

10.  Most  small-perturbation  theories  predict  ranges  of  flow  and  transport 
relations  under  which  the  bed-wave  amplitude  diminishes  with  time.  This  is 
interpreted  as  a  flat  bed. 

1 1 .  Riverine  bed  waves  produced  over  a  fairly  wide  range  of  wavelengths,  and 
the  resulting  bed  forms  are  the  superposition  of  waves  of  many  lengths 
(Jain  1971,  Jain  and  Kennedy  1974,  Nakagawa  and  Tsujimoto  1984). 
Free-surface  flows  are  observed  to  produce  two  spectral  peaks  in  the  bed 
waves  that  first  develop  on  the  initially  flattened  bed;  one  of  these  is  at  the 
wavelength  of  the  Airy  wave,  and  corresponds  to  the  bed  forms  produced 
by  standing  water-surface  waves.  Its  wavelength  is  given  by  the  Airy 
relation, 

U2  =  g/ktanhkd  (III.l) 

The  second  spectral-peak  is  produced  by  the  interaction  between  the  primary 
shear  flow  and  the  bed.  The  two  sets  of  bed-waves  interact,  by  the  shorter, 
faster-moving  ones  ovenaking  and  merging  with  the  slower,  longer  ones. 
This  "variance-cascade"  process  was  described  by  Jain  and  Kennedy  (1974) 
(see  Section  II.D.4  for  quotation  of  their  explanation). 

12.  The  spectral  peak  corresponding  to  the  Airy  celerity  is  not  present  in  bed 
forms  produced  by  closed-conduit  flows,  and  evolving  bed  forms  in  closed- 
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conduit  flows  do  not  exhibit  the  progressive  lengthening  described  by 
Raichlen  and  Kennedy  (1965)  and  analyzed  by  Jain  and  Kennedy  (1974). 
This  has  been  demonstrated  by  the  experiments  of  Nakagawa  and  Tsujimoto 
(1984). 

13.  Fully  developed  dune  spectra  exhibit  an  "equilibrium  range"  in  which  the 
spectral  density  function,  O,  varies  as 

(ni.2) 

where  k  =  2ji/L;  and  L  =  wavelength  (Hino  1968,  Jain  and  Kennedy  1974). 

For  ripple  spectra,  Hino  (1968)  suggests 

O-k-2  (ni.3) 

The  first  of  these  relations  has  been  well  corroborated  experimentally. 

1 4.  Some  data  compilations  indicate  that  very  long  bed  forms-megadunes— are 
a  distinctive  class  of  bed  waves.  This  is  suggested  by,  for  example,  the 
grouping  of  the  cluster  of  points  for  megadunes  in  the  lower-left  comer  of 
Fig.  B.4.3.  However,  spectral  density  functions  for  these  bed  forms  do  not 
indicate  that  they  are  statistically  different  from  other  bed  forms. 

15.  Megadunes  have  maximum  wavelengths  of  about  six-  to  eight-times  the 
flow  depth  (Jackson,  1976),  This  is  consistent  with  the  wavelength  of  the 
lower  end  of  the  equilibrium-spectrum  (O  -  k‘3)  range  (see  Section  II.D.2, 
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Eqs.  7  and  8),  as  has  been  shown  by  Hino  (1968).  This  corresponds  to  a 
vertical  line  in  Fig.  B.4.3,  which  is  seen  to  be  a  crude  estimate  of  dune 
wavelength. 

16.  A  better  approximation  to  the  wavelengths  of  dunes  can  be  obtained  from 
the  spectral  analyses  of  Nordin  (1971;  cited  in  Sec.  n.D.l)  and  Jain  and 
Kennedy  (1974;  see  Sec.  II.D.4)  indicating  that  better  estimates  of  the 
wavelength  are 


F2  = 


0.19 

kd 


(ni.4) 


((18)  in  Sec.  II.D.l);  and 


p5/3-Mil4 
2jc  ^ 


(ni.5) 


((19)  in  Sec.  II.D.4),  respectively. 

17.  Raudkivi's  (1963)  detailed  subcritical-flow  measurements  of  flow  over  a 
mobile-sand  and  a  rigid,  two-dimensional  model  of  a  ripple  (actually,  likely 
a  dune),  15-inches  long  and  1.1-inches  high,  revealed  several  interesting 
aspects  of  this  flow,  as  follow  (refer  to  Fig.  in.B.l): 

a.  Entrainment  and  bed-load  transport  take  place  on  the  stoss  slopes  of 
ripples  and  dunes  at  values  of  the  temporal  mean  bed  stress  well 
below  Shields'  critical  stress,  which  is  reached  more  than  half-way 
up  the  stoss  slope  of  the  bed  fonn. 
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b .  The  temporal  mean  shear  stress  is  a  maximum  on  the  crest  of  the  bed 
form,  and  there  is  very  nearly  equal  to  the  mean  shear  stress  for  the 
same  flow  over  same  sand  bed  in  its  initial  flat  condition.  This  flat¬ 
bed  value  of  shear  stress  appears  to  govern  the  amount  of  sediment 
transport. 

c.  The  streamlines  of  the  main  flow  are  approximately  sinusoidal, 
rising  over  the  wakes  and  falling  over  the  crests  of  ripples. 

d.  The  flow  over  the  wake  in  the  lee  of  a  ripple  reattaches  to  the  bed 
about  flve  to  eight  (average  of  about  six)  ripple-heights  downstream 
of  the  ripple  crest  At  the  reattachment  point,  the  bed  shear  stress  is 
practically  zero.  These  observations  are  consistent  with  those  made 
in  boundary  layers  downstream  from  impervious  "fences" 
perpendicular  to  walls;  and  downstream  from  negative  steps  on 
conduit  boundaries. 

e.  In  the  region  where  boundary  shear  stress  is  less  than  the  critical 
value,  the  sediment  entrainment  and  transport  are  due  primarily  to 
the  agitation  (intense  turbulence)  in  the  wakes  of  the  bed  forms. 

18.  There  is  presently  no  completely  reliable  predictor  for  the  conditions  of 
occurrence  of  the  different  bed  configurations  (ripples,  dunes/megadunes, 
flat  bed,  antidunes).  The  principal  stumbling  block  is  an  inadequate 
understanding  and  formulation  of  the  transport-rate  phase  shift  (see  Fig. 
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PRESSURE  IN  mm  OF  WATER 


Esuifis 


no.  4.-RIPPLE  FORMATION.  WATER  SURFACE  PROFILE,  AND  DISTRIBUTION  OF  PRESSURE 
DEVL\TIONS  AND  BED  SHEAR  STRESS 
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STRESS  To-loVf* 


FIG.  1.— Bed-Form  Chart  for  =  0.11  and  2.7-4.2  ld„  =  0.011  mm  and  0.088  mm) 


FIG.  2.~<6ed-Form  Chart  for  >  4.S-10  (d„  •  0.12  mmoO.200  mm) 


RG.  3.— Bed-Form  Chart  for  R,  >  10-16  (d„  =  0.15  iran-0.32  mm) 


III.B.2 
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FIG.  4.— Bad-Form  Chart  for  <■  IB-26  (d^  =  0.228  min-0.4S  mm| 


FIG.  5.— Bad-Form  Chart  for  B,  -  24-46  (d.  >  0.4  mm-0.57  mm) 


Figure  III.B.2.  (conrinued) 


157 


mk. 


TABLE  1. — Ranges  of  and  d„  for  Data  PtoRad  on  each  of  Bed>Form  Charts  Figs. 
2-7 


Figure 

number 

(1) 

Ranges  of  R. 

=  gti„)/v 

(2) 

Ranges  of  median  sediment 
size,  in  millimeters 

(3) 

1 

0. 1 1  and  2.7-4.2 

0.011  and  0.088-0.15 

2 

4.3-10 

0.12-0.20 

3 

10-16 

0.15-0.32 

4 

16-26 

0.23-0.45 

5 

24-50 

0.4-0.6 

6 

82-92  and  126-197 

0.93.  1.20.  1.35 

Table  m.l. 
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C.  Why  is  the  "Dune  Problem”  so  Intractable?  The  foregoing  outline  of  "what  is 
known",  and  the  selected-literature  survey  of  the  preceeding  chapter,  suggest,  accurately, 
that  the  mechanisms  responsible  for  the  inidation  and  small-amplitude  behavior  (the 
linearized  problem)  of  riverine  bed  forms  are  quite  well  understood  and  fairly  well 
formulated.  The  principal  remaining  stumbling  blocks  involve  the  usual  problems  attendant 
to  the  analysis  of  strongly  nonuniform  turbulent  shear  flows,  and  further  elucidation  and 
formulation  of  the  phase  shift  between  local  bed  displacement  and  local  sediment-transport 
rate.  The  perceptive  reader  will  have  noted,  however,  that  relatively  little  has  been  reported 
herein  concerning  analytical  description  of  the  properties  and  behavior  of  fully  developed 
bed  forms.  That  is  because  there  is  not  much  to  report,  either  analytical  or  experimental,  on 
this  problem.  The  notable  exceptions  are: 

1.  Fredsoe's  (1982)  analysis  of  the  shape  of  fully  developed  ripples  and 
dunes)  (see  Section  II.B.IO). 

2.  Haque  and  Mahmood's  (1985)  series  of  publications  reporting  analytical 
treatment  and  finite-element  simulation  of  flow  past  mature  ripples  and 
dunes  (see  Section  II.B.ll). 

3.  Several  dimensional-analysis  treatments  of  the  geometry  of  mature  bed 
forms  and  Vanoni's  (1974)  bed-form  charts. 

4.  Some  success  in  determining  "universal"  nondimensional  spectral  density 
functions  of  ripples  and  dunes  (see  Chapter  11,  Section  D). 

5.  Kennedy's  (1963)  finding  that  antidune  height  is  limited  by  the  stability  of 
the  accompanying  surface  waves  (sec  Section  n.B.3). 
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6.  The  relationship  between  the  geometiy  and  kinematics  of  ripples  and  dunes, 
and  sediment  discharge,  as  analyzed  by  Simons,  Richardson  and  Nordin 
(1965),  Kennedy  (1969),  and  Willis  and  Kennedy  (1978). 

In  an  attempt  to  answer  the  question  posed  in  the  heading  of  this  section,  let  us 
consider  what  a  complete  analysis  of  finite-amplitude  ripples  and  dunes  must  include,  and 
the  stumbling  blocks  encountered  in  its  development  We  will  consider  periodic,  fully 
developed  ripples  or  dunes  moving  with  celerity  c  without  changing  their  form,  as  shown 
in  Fig,  III.Cl.  These  waves  are  described  by 

T](x,t)  =  h  f  (x  -  ct)  (III.C.  1 ) 

where,  in  addition  to  the  nomenclature  defined  in  Fig.  III.C.  1,  f  =  function  of  (x-ct) 
defined  by  the  bed  profile  and  Eq.  III.C.  1,  and  t  =  time.  The  equation  of  sediment 
continuity  is 


at  at 


=  0 


(III.C.2) 


where  qs  =  volumetric  sediment  discharge  (which  includes  the  porosity).  Substitution  of 
Eq.  III.C.  1  into  Eq.  ni.C.2,  and  integrating  the  resulting  equation  over  x  yields 

qs(x,t)  =  cn(x,t)  +  C(t)  (III.C.3) 

where  C  =  constant  or  function  of  time,  and  represents  sediment  discharge  that  does  not 
participate  in  the  bed-form-migration  process. 
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In  Eq.  III.C.3,  "n  is  seen  to  vary  linearly  with  qs.  Thus  for  equilibrium  bed  forms, 
T|  must  evolve  such  that  it  produces  a  flow  field,  and  associated  turbulence  and  shear  field, 
that  at  every  point  along  the  bed  produces  a  q^  that  satisfies  Eq.  III.C.3.  This  relation  must 
be  satisfied  in  the  strongly  separated  flow  in  the  wakes  of  the  bed  forms;  in  the 
reattachment  zone,  where  the  separation  streamline  from  a  bed-form  crest  approaches  the 
stoss  slope  of  the  adjacent  downstream  bed  form;  and  in  the  nonuniform  internal  boundary 
layer  along  the  stoss  slope  downstream  of  the  reattachment  point.  In  short,  qs  must  be 
proportional  to  t]  throughout  the  complex  flow  above  ripples  and  dunes.  The  transport 
relation  must  also  take  into  account  the  effects  of  bed-form  slope  on  transport  rate, 
including  the  gravity  flow  of  sediment  down  the  lee  slopes.  It  must  also  consider  the 
strong  interchange  between  bed-load  and  suspended-load  discharges  that  occurs  along  the 
bed  profile,  including  the  wake  zone.  The  spectral  properties  of  the  bed  profile  would  have 
to  be  included,  especially  the  "production  of  bed  profile"  at  the  larger  wave  numbers,  and 
the  "decay  of  bed  profile"  at  smaller  wave  numbers.  Calculation  of  Tl(x,t)  and  the 
corresponding  flow  and  stress  fields  and  resulting  qs(x)  that  satisfy  these  relations,  as 
calculation  of  equilibrium  bed  forms  would  require,  is  indeed  a  tall  order.  To  appreciate  the 
difficulty  of  the  problem,  it  should  be  recalled  that  the  problem  of  calculating  the  total 
sediment  discharge  of  streams,  even  over  fiat  beds,  cannot  be  considered  as  satisfactorily 
solved. 
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IV.  A  NEW  THEORY  FOR  DUNE  HEIGHTS. 


This  model  is  concerned  primarily  with  megadunes,  of  the  type  described  by  Cany 
and  Keller  (1957),  among  others.  It  is  these  bed  forms  that  pose  the  greatest  problem  to 
navigation,  water-diversion  structures,  etc. 

The  model  is  based  on  the  following  ccmcepts  and  assumptions: 

1 .  The  grain-roughness  shear  stress,  t’,  varies  linearly  along  the  stoss  face  of 
each  dune  (see  Fig.  m.C.l),  from  zero  a  distance  oh  downstream  from  the 
dune  crest  to  the  flat-bed  value  (based  on  the  mean  depth  d)  To>  dune 
crest.  This  assumption  is  well  supported  by  the  detailed  measurements  of 
Raudkivi  (1963)  of  velocity  and  stress  distributions  in  flows  over  dunes. 
Based  on  his  experiments,  and  the  work  of  others  cited  by  him,  a  s  5.  The 
negative  t'  upstream  from  this  point  to  the  dune  crest  will  be  neglected. 
Thus,  the  average  energy  slope  due  to  the  grain  roughness.  S’,  is  expressed 
by 


pgdS'  =^Xxo 


L-ah 

L 


or 


S’=  1.2|^A.foF2(l-a^)  (IV.  1) 

where,  in  addition  to  the  terms  defined  in  Fig.  III.C.l,  X  =  coefficient  (=  1); 
fo  =  rigid-flat-bed  Darcy- Weisbach  friction  factor,  and  F  =  Froude  number. 
The  factor  1.2  is  introduced  to  take  into  account  Karim  and  Kennedy's 
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(1989)  finding  that  the  friction  factors  of  flows  over  mobile  flat  beds  are 
about  20  percent  greater  than  those  of  flows  over  rigid  beds  with  equal 
roughness  size. 


2.  The  energy  slope  produced  by  form  drag  on  the  dunes  will  be  estimated 
from  Carnot's  formula  for  the  head  loss  across  an  abrupt  expansion  in  a 
conduit,  as  was  done  by  Engelund  (1966).  The  average  energy  slope,  S", 
produced  by  this  expansion  loss  is  (refer  to  Fig.  IV.  1) 


where  f  =  Darcy-Weisbach  friction  factor  and  F  =  Froude  number,  (rv.3) 
yields 


av.4) 


4.  The  ratio  ^  will  be  estimated  from  the  results  of  the  spectral  analyses  of 

Nordin  (1971)  and  Jain  and  Kennedy  (1974)  (see  item  16  of  Sec.  III.B). 
Specifically, 


F2  =  Ci_QL 

dk-2;,d 


(IV.5) 


will  be  adopted.  It  is  expected  that  Ci  s  0.25  for  megadunes.  (VI.S)  is  a 
reasonable  approximation  to  the  "megadune"  points  in  the  lower  left  of  Hg. 
B.3.2. 


Karim  and  Kennedy's  (1989;  their  Eq.  17)  regression  relation  will  be 
utilized  for  the  friction-factor  ratio  f/fo: 


U 

'g(s-l)D5o 


=  6.683  (^)  -626  (S)  .503  (^)  -.465  (IV.6) 


where  s  =  specific  gravity  of  bed  material;  and  D50  =  median  bed-particle 


6.  The  rigid-flat-bed  friction  factor,  fo,  is  taken  from  Karim  and  Kennedy 
(1989)  to  be 
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6.25  +  2.5in 


2.5  D50J 


(IV.7) 


7.  (IV.4)  and  (IV.5)  yield 


f  .  . 


where  fo  is  given  by  (IV.7);  Ci  ~  0.25;  X  r  1;  Cd  ~  1;  and  7-  is  given  by 
(IV.6). 


.  l.2Xafo  ,  1.2 


8.  (IV.8)  can  be  simplified  (by  neglecting  the  terms  in  -g — ^ 

which  are  small  for  all  but  flows  over  nearly  flat  beds)  to 


and  X , 


h  F-v/  2jrf 
d”2  VCdCi 


(IV.9) 


Several  comments  concerning  the  prospective  usefulness  of  (rV.9)  and  (IV.  10)  are 


in  order 


They  valid  only  for  flows  that  produce  ripples  and/or  dunes.  Karim  and 
Kennedy  (1989)  suggest  that  this  flow  regime  occurs  over  the  range 


0.06  < - - - <1.3 

pg(s-l)  D50 


(IV.  10) 


where  0.06  corresponds  to  the  Shields  dimensionless  shear  stress  for 
incipient  motion.  The  limits  given  in  (FV.IO)  are  tentative  and  approximate. 
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Additionally,  dunes  can  be  expected  only  if  the  Froude  number  is  less  than 
about  0.6  to  0.7.  Vanoni’s  (1974)  bed-form  charts  (see  item  18  of  Sec. 
ni.B,  and  Fig.  III.B.2)  also  should  be  consulted  to  determine  if  dunes  can 
be  expected. 

2.  Because  ripples  and  dunes  are  composed  of  spectra  of  wavelengths  and 
amplitudes,  any  prediction  for  these  quantitides,  or  their  ratio,  should  be 
regarded  as  yielding  representative  values. 

3 .  (rv.9)  can  be  expected  to  be  valid  only  for  fairly  large  values  of  ^  . 
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V.  VERIFICATION 


The  relative  dune-height,  predictions  of  six  of  the  empirically  based  relations 

were  compared  with  data  reported  by  Guy  et  al.  (1966).  Only  their  runs  conducted  in  the 
eight-foot-wide  flume  which  had  dunes  as  the  reported  bed  configuration  were  used  in  this 
comparison.  Regretably,  no  field  data  collection  suitable  for  verification  were  found. 

The  ranges  of  variables  covered  by  these  experiments  utilized  in  the  comparison  are: 

d  =  flow  depth  (ft):  0.30  -  1.33 
O50  =  grain  size  (mm):  0. 19  -  0.93 
F  =  Froude  number  0.25  -  0.65 
S  =  Slope:  0.37  x  10-3 . 4  37  x  10-3 
V  =  velocity  (ft/s):  1:30  -  3.32 

The  following  relations  were  used  in  the  verifications: 

(1)  Yalin  (1964) 


h 

d" 


'k  _  Qc(ps  -  P)gD 
T  pgSd 


AD0C 

Sd 


(9c  =  0.06) 


169 


(2)  GiU(1971) 


h  1  -  F2  /  V 
^  2na 


(assume  f  is  independent  of  U) 


a  =  0.5  (triangular  dunes) 


n  =  3/2  (corresponds  to  Meyer-Peter  and  Mueller  equation). 


(3)  Ranga  Raju  and  Soni  (1976) 


^=6.5xI03-§^(t;)S/3 

FjF2 


Fl=7%r-  F2 

VgR^ 


'  Ap 
—  gD 

P 


ns=^:^(m 


Uinm/s; 


(4)  Fredsoc  (1982) 


de  de 
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Fig.  B.104  portrays  this  relation  after  0,  <l>b,  <I>s  are  expressed  as  functions 
of0: 


0c  =  0.06 

The  results  of  this  verification  are  presented  in  Figs.  V.l  to  V.5.  These  figures 
suggest  that  none  of  the  relations  evaluated  is  fully  satisfactory.  The  model  developed 
herein  appears  to  give  the  best  results,  and  could  be  improved  by  calibration  to  adjust  the 
constants  it  contains. 

(6)  Present  paper  (Chapter  FV). 

h  l/l.2Xafo  rn.2XafoV  27tF2fo  ^f  1.2X>|]i/2l 

d "  2 1  8Cd  ^  ^  CdCi  Uo  ■  2  jJ  I 


where 


8 


|6.25  +  2.5  Jtn 


a  =  5 
X=1 


Cl  =  0.25 
Cd=1 


d...  f 

2.5  D50J 
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MEASURED  RELATIVE  DUNE  HEIGHTS,  h/d 


MEASURED  RELATIVE  DUNE  HEIGHTS,  h/d 


MEASURED  RELATIVE  DUNE  HEIGHTS,  h/d 


Figure  V.3.  Verification  of  the  dune-height  prediction  of  Ranga-Raju  and  Soni  (1971). 
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MEASURED  RELATIVE  DUNE  HEIGHTS,  h/d 


6 


MEASURED  RELATIVE  DUNE  HEIGHTS,  h/d 


MEASURED  RELATIVE  DUNE  HEIGHTS,  h/d 


VI.  DUNE-CONTROL  METHODS 


In  the  extensive  literature  search  carried  out  in  the  course  of  this  study,  no 
innovative  means  of  dune  control  were  found.  Dredging  to  obtain  shon-term  relief  of  local 
problems  seems  presently  to  be  the  common  practice. 

A  technique  whose  potentials  as  a  control  measure  have  not  been  fully  explored  is 
the  submerged  vanes  technique.  Both  laboratory  and  field  tests  (Odgaard  and  Kennedy 
1983,  Odgaard  and  Spoljaric  1986,  Odgaard  and  Mosconi  1987,  Wang  1989,  Fukuoka 
1989,  Fukuoka  and  Watanabe  1989)  suggest  that  this  technique  has  a  broad  range  of 
applications  including  dune  control.  The  vanes  are  small  flow-training  structures  (foils) 
designed  to  modify  the  near-bed  flow  pattern  and  redistribute  flow  and  sediment  transport 
within  the  channel  cross  section.  The  structures  are  installed  at  an  angle  of  attack  of  15-25° 
with  the  flow,  and  their  initial  height  is  0.2-0.4  times  local  water  depth  at  design  stage. 
The  vanes  function  by  generating  secondary  circulation  in  the  flow.  The  circulation  alters 
magnitude  and  direction  of  the  bed  shear  stress  and  causes  a  reduction  in  velocity  and 
sediment  transpon  in  the  vane  controlled  area.  As  a  result,  the  bed  configurations  change. 
The  river  bed  aggrades  in  the  vane  controlled  area  and  degrades  outside. 

The  technique  was  originally  developed  to  stop  or  reduce  bank  erosion  in  river 
curves.  In  this  application  the  vanes  are  laid  out  so  that  the  vane  generated  secondary 
current  eliminates  the  centrifugally  induced  secondary  current,  which  is  the  root  cause  of 
bank  undermining.  The  first  field  test  with  this  application  was  in  a  bend  of  East 
Nishnabotna  River,  Iowa  (Odgaard  and  Mosconi  1987).  The  results  from  this  test  are  very 
satisfactory.  Other  successful  field  tests  are  reported  by  Fukuoka  and  Watanabe  (1989) 
and  Kunzig  (1989). 
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Recent  laboratory  and  field  tests  suggest  that  the  vane  technique  is  also  effective  in 
ameliorating  shoaling  problems  in  rivers  (Odgaard  and  Wang  1990).  Figures  VI.  1  and 
VI.2  show  results  of  a  laboratory  experiment  in  which  vanes  were  installed  along  the  right 
side  of  a  straight  laboratory  channel.  The  objective  was  to  generate  an  increase  in  flow 
depth  in  the  left  portion  of  the  channel.  The  vanes  were  0.8  mm  thick  sheet  metal  plates, 
7.4  cm  high  and  15.2  cm  long,  installed  in  arrays  with  four  vanes  in  each  array  and  angled 
20°  toward  the  bank.  Average  flow  depth  was  about  17  cm  and  velocities  were  of  the 
order  of  0.4  m/s.  It  is  evident  that  this  vane  system  caused  a  considerable  redistribution  of 
sediment.  The  vanes  reduced  the  depth  near  the  right  bank  by  about  50%.  This  caused  the 
depth  near  the  left  bank  to  increase  by  20-30%. 

One  of  the  most  important  observations  made  in  these  tests  was  that  the  vane 
induced  changes  in  cross-sectional  profile  occurred  without  causing  measurable  changes  of 
the  area  of  the  cross  sections  and  of  the  longitudinal  slope  of  the  water  surface.  This 
observation  is  important  because  it  implies  that  the  vanes  will  not  cause  any  changes  of  the 
stream's  sediment  transpon  capacity  upstream  or  downstream  from  the  vane  field  and 
therefore  should  not  alter  the  overall  characteristics  of  the  stream. 

Another  important  observation  was  that  the  bedforms  in  the  deepened  part  of  the 
channel  did  not  appear  to  be  larger  than  before  vanes  were  installed.  It  appeared  that  the 
vanes  lowered  the  bed  level  outside  the  vane  field  without  causing  the  height  of  the 
bedforms  there  to  be  increased.  This  observation  suggests  that  vanes  may  be  an  effective 
means  of  controlling  not  only  the  depth  but  also  the  size  of  the  bedforms.  However,  more 
experimentation  is  needed  to  quantify  the  performance  of  the  vanes  in  this  regard. 

It  is  recommended  that  the  use  of  submerged  vanes  for  control  of  bedforms  be 
investigated  further.  A  logical  approach  would  be  to  use  the  theoretical  basis  obtained 
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herein  (the  relationships  between  bedform  geometry  and  flow  and  sediment  characteristics) 
for  designing  and  conducting  a  series  of  experiments  in  a  straight,  movable-bed  hydraulic 
model.  The  main  objective  would  be  to  develop  relationships  between  vane-system 
characteristics  (design  and  layout),  flow  and  sediment  conditions,  and  dune  dimensions 
within  and  outside  the  vane  controlled  region  of  the  model.  The  experiments  should  be 
designed  so  that  the  feasibility  of  the  technique  in  typical  prototype  simations  can  be 
determined.  The  model  could,  for  example,  be  a  representation  of  a  prototype  problem 
reach.  Other  techniques  may  also  emerge  from  such  experiments. 
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Figure  VI.  I .  Submerged  vanes  for  sediment  management:  test  in 
straight,  alluvial  channel.  Upstream  view  of  nearly 
drained  channel. 
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DEPTH  BELOW  WATER  SURFACE,  d/d 


1.5 


Figure  VI.2.  Submerged  vanes  for  sediment  management:  test  in 

straight  alluvial  channel.  Velocity  and  depth  distributions 
with  and  without  vanes. 
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VII.  SUMMARY  AND  CONCLUSIONS 


The  principal  findings  forthcoming  from  the  literature  survey  are  presented  in 
Section  III.B,  which  should  be  considered  a  part  of  the  summary  and  conclusions  of  this 
study. 


Five  of  the  published,  empirically  based  predictors  for  relative  dune  height,  h/d, 
were  evaluated  (Chapter  V).  None  of  them  was  found  to  be  satisfactory.  A  new  dune- 
height  predictor  was  derived  (Chapter  IV),  using  in  "inverse  approach",  in  which  the  dune 
height  was  estimated  from  a  verified  predictor  for  alluvial-river  friction  factors.  This 
relation  was  found  to  give  more  accurate  predictions  of  h/d  than  the  previously  published 
ones.  Its  accuracy  no  doubt  could  be  improved  by  calibration  against  other  data  to  refine 
the  constants  it  contains;  and  by  omitting  some  of  the  assumptions  made  in  its  formulation, 
although  this  would  be  done  at  the  expense  of  greater  computational  complexity. 

The  most  promising,  innovative  means  of  dune  control  is  judged  to  be  use  of 
submerged  vanes,  as  described  by  Odgaard  and  Wang  (1990).  As  a  matter  of  fact,  no 
other  innovative  means  was  discovered.  There  is  a  growing  body  of  field  experience 
which  demonstrates  that  submerged  vanes,  suitably  deployed,  are  very  effective  sediment 
management  and  control  devices.  It  is  recommended  that  their  use  for  control  of 
objectionable  shoaling  resulting  from  dunes  be  further  explored. 
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